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Abstract 

Extending investigations of Barker, Humpherys, Lafitte, Rudd, and Zumbrim for com- 
pressible gas dynamics and Freistiihler and Trakhinin for compressible magnetohydrody- 
namics, we study by a combination of asymptotic ODE estimates and numerical Evans 
function computations the one-dimensional stability of parallel isentropic magnetohydrody- 
namic shock layers over the full range of physical parameters (shock amplitude, strength 
of imposed magnetic field, viscosity, magnetic permeability, and electrical resistivity) for a 
7-law gas with 7 g [1,3]. Other 7- values may be treated similarly, but were not checked 
numerically. Depending on magnetic field strength, these shocks may be of fast Lax, in- 
termediate (overcompressive), or slow Lax type; however, the shock layer is independent 
of magnetic field, consisting of a purely gas-dynamical profile. In each case, our results 
indicate stability. Interesting features of the analysis are the need to renormalize the Evans 
function in order to pass continuously across parameter values where the shock changes type 
or toward the large-amplitude limit at frequency A = and the systematic use of winding 
number computations on Riemann surfaces. 

1 Introduction 

In this paper, continuing investigations of [31 127^ [T3]. we study by a combination of asymptotic 
ODE estimates and numerical Evans function computations the one-dimensional stability of 
parallel isentropic magnetohydrodynamic (MHD) shock layers over a full range of physical 
parameters, including arbitrarily large shock amplitude and strength of imposed magnetic field, 
for a 7-law gas with 7 G [1,3], with our main emphasis on the case of an ideal monatomic or 
diatomic gas. The restriction to 7 G [1,3] is an arbitrary one coming from the choice of 
parameters on which the numerical study is carried out; stability for other 7 can be easily 
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checked as well. (Note that our analytical results are for any 7 > 1.) In each case, we obtain 
results indicative of stability. Recall that Evans stability, defined in terms of the Evans function 
associated with the linearized operator about the wave, by the "Lyapunov-type" results of 
[381 l39l [5H [52I [23l [Ml [36] , implies linear and nonlinear stability for all except the measure-zero 
set of parameters on which the characteristic speeds of the endstates coincide with the shock 
speed or each other0 

Parallel shocks may be of fast Lax, intermediate (overcompressive) , or slow Lax type de- 
pending on magnetic field strength; however, the shock layer is independent of magnetic field, 
consisting of a purely gas-dynamical profile. Thus, the study of their stability is both a natural 
next step to and an interesting generalization of the investigations of stability of gas-dynamical 
shocks in |27j . See also the investigations of stability of fast parallel Lax shocks in certain param- 
eter regimes in [TSj using energy methods, and of general fast Lax shocks in the small-magnetic 
field limit in |20[ [T9] using Evans function techniques. 



1.1 Equations 

In Lagrangian coordinates, the equations for compressible isentropic MHD in one dimension 
take the form 

Vt - Ulx = 0, 

Ult + {p+ {1/2ixq){BI + Bl))x = (((2/i + r])/v) 

U2t - {{l/^io)B\B2)x = {{lJ'/v)u2x)x, 

U3t - {{l/flo)BlBs)x = {ifx/v)u3x)x, ^ ■ ^ 

{vB2)t - {Blu2)x = i{l/(^tJ'0v)B2x)x, 

, {vB3)t - {Blu3)x = {{l/criJ-ov)B3x)x, 

where v denotes specific volume, u = (ui,?i2,U3) velocity, p = p{v) pressure, B = {Bl, B2, B3) 
magnetic induction, B^ constant, and /i > and r] > the two coefficients of viscosity, /io > 
the magnetic permeability, and a > the electrical resistivity; see [3l [lOl [331 ES] for further 
discussion. 

We restrict to an ideal isentropic polytropic gas, in which case the pressure function takes 
form 

p(v) = av-^ (1.2) 

where a > and 7 > 1 are constants that characterize the gas. In our numerical investigations, 
we shall focus mainly on the most common cases of a monatomic gas, 7 = 5/3, and a diatomic 
gas, 7 = 7/5; more generally, we investigate all 7 G [1,3]. With brief exceptions (e.g.. Section 
[13]), we take 

T] = -2fi/3, (1.3) 

as typically prescribed for (nonmagnetic) gas dynamics [5]. 

Here, we are allowing u and B to vary in full three-dimensional space, but restricting spatial 
dependence to a single direction ei measured by x. That is, we consider planar solutions, 

^For these degenerate cases, the stabihty analysis has not been carried out in the generahty considered here. 
However, see the related analyses for Lax shock of [251122] in the case that shock and characteristic speed coincide 
and [51] in the case that characteristic speeds coincide, which suggest that the shocks may be nonetheless stable. 
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or three-dimensional solutions with one-dimensional dependence on spatial variables. Note 
that the divergence- free condition div^B = of full MHD reduces in the planar case to our 
assumption that Bi = constant = B*. In the simplest, parallel case 

B2 = B3 = 0; U2=us = 0, (1.4) 

equations (jl.lh reduce to the one-dimensional isentropic compressible Navier-Stokes equations 

vt - = 0, . 
uit+Px = {{{2^x + r|)/v)ulx)x■ 
l\l the remainder of the paper, we study traveling-wave solutions in this special parallel case 
and their stability with respect to general (not necessarily parallel) planar perturbations. 



1.2 Viscous shock profiles 

A viscous shock profile of (jl.ip is an asymptotically-constant traveling-wave solution 
{v,u, B){x,t) = {v,u, B){x — st), lim =iv±,u±,B±). 

In the parallel case, these are of the simple form 

{v,u,B){x - St) = {v,ui,0,0,Bl,0,0){x - st), 
where {v,ui) is a gas-dynamical shock profile satisfying the traveling- wave ODE 

-SVx - Ulx = 0, 

-suix +Px = (((2/i + r])/v)uix)x- 



(1.6) 



(1.7) 



1.3 Rescaled equations 

By a preliminary rescaling in x, t, we may arrange without loss of generality = 1. Following 
the approach of [271 ISHl |28] , we now rescale 

{v,ui,U2,U3,fio,x,t,B) -, , — , — ,e/io, -es{x - st),es t, — 

\e es e e s 

holding ^, cr fixed, where £ := V-, transforming (jl.ip to the form 

Vt+Vx- Uix = 

1 



Ult +uix+ [ av ' + 



[Bl + Bl] 



1 



"3t + - i—BlBy, 

{VB2)^ + {VB2), ~ {BIU2) 

ivBs)^ + (vB,)^ - iB*us) 



U2x 
V 



\ V Ji 



1 



1 



B2x 

Bzx 



(1.8) 
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where p{v) = aov~^ and a = aoe~'^~^ s'"^ . 

By this step, we reduce without loss of generahty to the case of a shock profile with speed 
s = — 1, left endstate 

{v,uuU2,U3,Bi,B2,Bs)^ = (1,0, 0,0, Si*, 0,0), (1.9) 

and right endstate 

iv,ui,U2,U3,Bi,B2,Bs)+ = iv+,v+ - 1,0, 0,^1*, 0,0), (1.10) 

satisfying the profile ODE 

{2fi + 7]^ = H{v, v+) := v{v-l + aiv-"^ - 1)) (1.11) 

(obtained by integrating (jl.7p and substituting the first equation into the second) where 1 = 
V- > > and (setting v' = at v = v+ and solving) 

a = --^ = .,^. (1.12) 

See [H EZj for further details. 

Proposition 1.1 ([4J). For each 7 > 1, < t;+ < 1 — e, e > 0, (jl.lip has a unique (up to 
translation) monotone decreasing solution v decaying to its endstates with a uniform exponential 
rate, independent of Vj^, 7. In particular, for < f + < and v{Q) := v+ + j^, 

\v{x) -v+\ < (^-^^e'^ x>0, (1.13a) 
\v{x) -V-\ < i^-je— x<0. (1.13b) 



Corollary 1.2. Initializing v{0) := v+ + as in Proposition \l.ll v converges uniformly as 



1— tanh I 



^ to a translate vq of \^(2f^ ^ 

Proof. By (|1.12p . a ~ ^ as f+ ^ 0, whence the result follows on any bounded set 
|x| < L by continuous dependence, taking the limit as a ^ in (jl.lip to obtain a limiting 
flow of v' = ^2;j+?j'* • Taking now L — > 00, the result follows for > M by t;+ ^ and 
\v-v+\< Ce-^^'; see ifTTSall - ffLTSb]) . □ 



1.4 Families of shock profiles 

At this point, we have reduced our study of parallel shock stability, for a fixed gas constant 7, 
to consideration of a one-parameter family of profiles indexed by the right endstate 1 > W4. > 
and a four-parameter family of equations (jl.Sp indexed (through (jl.l2p ) by Vj^ and the three 
remaining physical parameters 

/xo > 0, cj > 0, ^1 > 0, (1.14) 
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where we have taken without loss of generality to be nonnegative by use of the symmetry 
under B — > —B of (jl.ip . Here, the small-amplitude limit corresponds to — > u _ = 1 and the 
large-amplitude limit to ^ 0, where in this scaling the amplitude is given by |u_ — 

A straightforward computation shows that the characteristics of the first-order hyperbolic 
system obtained by neglecting second-derivative terms in p.8|) at the endstates v± have values 

where c{v) := y/ —p'{v) = \J ^av^^^^ is the gas-dynamical sound speed, satisfying c+ > 1 > c_. 
Thus, the shock is a Lax 1-shock for < Bl < ^^qv+, meaning that it has six positive 
characteristics at and one at f+; an intermediate doubly overcompressive shock for ^/iof+ < 
^1 < meaning that it has six positive characteristics at and three at and a Lax 

3-shock for ^/Jiq < B^, meaning that it has 4 positive characteristics at and three at v^. 

For Lax 1- and 3-shocks, the profile ()1.6p is generically (and always for 1-shocks) unique 
up to translation as a traveling-wave solution of the full equations connecting endstates (jl.9p 
and (jl.lOp . i.e., even among possibly nonparallel solutions. That is, it lies generically within 
a one-parameter family {f7^} = {{v,ui,ii2,U3, B2, Bs)^} of viscous shock profiles, ^ G M, with 
U((x) := U{x - 0- For overcompressive shocks, it lies generically within a three-parameter 
family {({), ui,U2,U3, B2, B^)^} of viscous profiles and their translates, ^ G M^, of which it is the 
unique parallel solution up to translation [38j. For further discussion of hyperbolic shock type 
and its relation to existence of viscous profiles, see, e.g., [37 1 l55 l \50 \ [STl [38] . 

1.5 Evans, spectral, and nonlinear stability 

Following [55l EHl E] , define spectral stability as nonexistence of nonstable eigenvalues 5RA > 
of the linearized operator about the wave, other than at A = (where there is always an eigen- 
value, due to translational invariance of the underlying equations). A slightly stronger condition 
is Evans stability, which for Lax or overcompressive shocks may be defined [551 [38l [28] as nonva- 
nishing for all 5RA > of the Evans function associated with the integrated eigenvalue equation 
about the wave. See [U [TBI EOl [SU [38] for a general definition of the Evans function associated 
with a system of ordinary differential equations; for a definition in the present context, see 
Section [2j Recall that zeros of the Evans function (either integrated or nonintegrated) agree 
with eigenvalues of the linearized operator about the wave on {KA > 0} \ {0}, so that Evans 
stability implies spectral stability. 

The following "Lyapunov-type" result of Raoofi |l5], specialized to our case, states that, for 
generic parameter values, Evans stability implies nonlinear orbital stability, regardless of the 
type of the shock; see also [39] [5T| [M] [i6] . 

Proposition 1.3 ([^). Let U := (v,ui,U2,U3, B2, B-^) be a parallel viscous shock profile of 
(|l.ip - ()1.2p connecting endstates ()1.9p - (|1.10p . with characteristics (jl.lSp distinct and nonzero, 
that is Evans stable. Then, for any solution U := {v,ui,U2,U3, B2, B^) of (II. ip with L^ n 
initial difference and L^ -first moment Eq := ||C/(-, 0) — [/||^ip|j:^3 and Ei := |||x| |C/(-, 0) — t/| ||^i 
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sufficiently small and some uniform C > 0, U exists for all t > 0, with 

\\U{;t) - U{- - st)\\LinH3 < CEo (stability). (1.16) 
Moreover, there exist a{t), Ooo such that 

\\if{.,t)-U^'^'\--stmLP <CEo{l + t)-^'/'^^'-'M, (1.17) 

and 

|a(t)-aoo|, {l + t)^/^\a{t)\ < C(e) max{^o, ^i}(l + t)"^/^+", (1-18) 
for all 1 < p < OO, e > arbitrary (phase-asymptotic orbital stability). 

Finally, recalling that Evans stability for Lax shocks is equivalent to the three conditions 
of spectral stability, transversality of the traveling wave as a connecting orbit of (jl.lip . and 
inviscid stability of the shock while Evans stability for overcompressive shocks is equivalent to 
spectral stability, transversality, and an "inviscid stability" -like low- frequency stability condi- 
tion generalizing the Lopatinski condition of the Lax case [551 138| [5T] , we obtain the following 
partial converse allowing us to make stability conclusions from spectral information alone. 

Proposition 1.4. A parallel viscous shock profile of (jl.ip - (jl.2p . (jl.9p - (jl.l0p . that is a Lax 
1-shock and spectrally stable is also Evans stable (hence, for generic parameters, nonlinearly 
orbitally stable). A parallel viscous shock profile that is an intermediate (overcompressive) 
shock, spectrally stable and low-frequency stable is Evans stable. A parallel viscous shock profile 
that is a Lax 3-shock, spectrally stable, and transverse is Evans stable. For /i = 1 and B\ > 

max I \J-^, \[i } ~'~ -^'^^ 3-shocks are transverse. For parallel viscous shocks of any type, 
spectral stability implies Evans (and nonlinear) stability on a generic set of parameters. 

Proof. Lax 1-shocks and intermediate-shocks, as extreme shocks (i.e., all characteristics entering 
the shock from the — oo side), are always transversal |38j. One-dimensional inviscid stability 
of either Lax 1- or 3-shocks follows by a straightforward calculation using decoupling of the 
linearized equations into {v,ui) and (^2,-62) and (u^jB^) systems |6lll9l[13]. Transversality 
for large is shown in Proposition IB. 31 Finally, both transversality (in the Lax 3-shock 
case) and (in the overcompressive case) low-frequency stability conditions can be expressed 
as nonvanishing of functions that are analytic in the model parameters, hence either vanish 
everywhere or on a measure zero set. It may be shown that these are both nonvanishing for 
sufficiently weak profiles |1 — v+| smallH hence they are generically nonvanishing. From these 
facts, the result follows. □ 



^ For Lax 3-shocks, transversality follows for small amplitudes by the center-manifold analysis of [32] ■ For 
overcompressive shocks, taking Bl = (l/2)(y^/Io + v^/J^T) «+ 1, using decoupling of and {uj,Bj) 

equations and performing a center manifold reduction in the {uj, Bj) equation of the traveling- wave ODE written 
as a first-order system, j — 2, 3, we find that this reduces in each case to a one-dimensional fiber, whence decaying 
solutions of the linearized profile equation, corresponding to variations other than translation in the family of 
profiles , are of one sign and thus have nonzero total integral J^°°{uj, Bj){x)dx. But this is readily seen |50| 
to be equivalent to low-frequency stability in the small-amplitude limit. 
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Remark 1.5. Our numerical results indicate Evans stability for all parameters, which implies 
in passing uniform transversality of 3- and overcompressive-shock profiles and low-frequency 
stability of overcompressive profiles. Transversality is a minimal condition for orbital stability, 
being needed even to guarantee existence of the smooth mamjold tj°' under discussion J3^ . As 
discussed above, it is not implied by spectral stability aloneo 



1.6 The reduced linearized eigenvalue equations 

Linearizing (II. Sp about a parallel shock profile {v, ui, 0, 0, 3^,0, 0), we obtain a decoupled system 

Vt + Vx- Ulx = 



U2t + U2x 



aj yv 
1 

1 



-7-1, 



1 



(2^ + 7? — + —V 



{BlB2)x 



ust + usx {BlB-i)x 



ivB2), + ivB2) 
ivBs)^ + (vB^) 



iBlu2)x 
{Btus)x 



U2x 

V 
U3x 

V 

1 

1 

CJ^O 



(1.19) 



B' 



2x 



V 

Bsx 

V 



consisting of the linearized isentropic gas dynamic equations in (v, ui) about profile (v, ui), and 
two copies of an equation in variables {uj,vBj), j = 2,3. 

Introducing integrated variables V := J v, U := f ui and Wj := / uj, aj := / vBj, j = 
2,3, we find that the integrated linearized eigenvalue equations decouple into the integrated 
linearized eigenvalue equations for gas dynamics in variables {V, U) and two copies of 



\w + w' 



Xa 



a 



Bla' 
B^w 



w 
i) 
1 



(1.20) 



in variables {wj, Uj), j = 2, 3. 

As noted in [55l [38l [2H] > spectral stability is unaffected by the change to integrated variables. 
Thus, spectral stability of parallel MHD shocks, decouples into the conditions of spectral stabil- 
ity of the associated gas-dynamical shock as a solution of the isentropic Navier-Stokes equations 
(II. sp . and spectral stability of system (ll.20p . Assuming stability of the gas-dynamical shock (as 
has been verified in great generality in [271 [29]), spectral stability of parallel MHD shocks thus 
reduces to the study of the reduced eigenvalue problem (ll.20p . into which the shock structure 



■^Thus, for example, the spectral stability results obtained by energy estimates in [13 for intermediate- or 
Lax 3-shocks do not by themselves imply linearized or nonlinear stability, but require an additional study of 
transversality /low- frequency stability. 
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enters only through density profile v. Likewise, the Evans function associated with the full 
system (|1.19p decouples into the product of the Evans function for the gas-dyamical eigenvalue 
equations and the Evans function for the reduced eigenvalue problem p.20p . Thus, assuming 
stability of the associated gas-dynamical shock, Evans stability of parallel MHD shocks reduces 
to Evans stability of (|1.2U|) . 

Remark 1.6. The change to integrated coordinates removes two additional zeros of the Evans 
function for the reduced equations (jl.20p that would otherwise occur at the origin in the over- 
compressive case, making possible a unified study across different parameter values/shock types. 



1.7 Analytical stability results 

1.7.1 The case of infinite resistivity /permeability 

We start with the observation that, by a straightforward energy estimate, parallel shocks are 
unconditionally stable in transverse modes (n, B) in the formal limit as either electrical resis- 
tivity a or magnetic permeability fiQ go to infinity, for quite general equations of state. This is 
suggestive, perhaps, of a general trend toward stability. 

Theorem 1.7. In the degenerate case fiQ = oo or a = oo, parallel MHD shocks are transversal, 
Lopatinski stable (resp. low- frequency stable), and spectrally stable with respect to transverse 
modes {u,B), for all physical parameter values, hence are Evans (and thus nonlinearly) stable 
whenever the associated gas-dynamical shock is Evans stable. 

Proof. By Proposition 11.41 Lopatinski stability holds for Lax-type shocks, and transversality 
holds for Lax 1-shocks and overcompressive shocks. Noting that the (decoupled) transverse 
part of linearized traveling- wave ODE for u = oo or /^q = oo reduces to {d — 1) copies of the 
same scalar equation, and recalling that transversality /Lopatinski stability hold always for the 
decoupled gas-dynamical part [27j . we readily verify low- frequency stability in the overcom- 
pressive case and transversality in the Lax 3-shock case as Thus, we need only verify 
transverse spectral stability, or nonexistence of decaying solutions of (ll.20p . 
For o" = oo, we may rewrite ()1.20p in symmetric form as 

liovXw + fiQVw' - Bla = fifiow", 
Aa -|- q' — Blw' = 0. 

Taking the real part of the complex L^-inner product of w against the first equation and a 
against the second equation and summing gives 



^ For scalar equations, transversality is immediate. Likewise, decaying solutions of the linearized profile 
equation, corresponding to variations other than translation in the family of profiles U" , are necessarily of one 
sign and thus have nonzero total integral J^°° {uj , Bj){x)dx. But this is readily seen ^D] to be equivalent to 
low-frequency stability in the small-amplitude limit. 
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a contradiction for ?R.X > and w not identically zero. If w = on the other hand, we have a 
constant-coefficient equation for a, which is therefore stable. The ^0 = 00 case goes similarly; 
see Appendix IB. 2[ □ 

Notably, this includes all three cases: fast Lax, over compressive, and slow Lax type shock. 
Further, the same proof yields the result for the more general class of equations of state p(-) 

p{v+)-p{v-) 

D-l- —V— 

obtain in particular the following asymptotic results. 



satisfying v^-v — < Oj so that Vx < for s < 0. With the analytical results of [27], we 



Corollary 1.8. For a = 00 or fiQ = 00, parallel isentropic MHD shocks with ideal gas equation 
of state, whether Lax or overcompressive type, are linearly and nonlinearly stable in the small- 
and large- amplitude limits f + ^ 1 and f+ — > 0, for all physical parameter values. 

1.7.2 Bounds on the unstable spectrum 

By a considerably more sophisticated energy estimate, we can bound the size of unstable eigen- 
values uniformly in 1 > > and the gas constant 7 > 1 to a ball of radius depending on a, 
/io, Bl, a crucial step in studying the limit f+ 0. 

Theorem 1.9. Nonstable eigenvalues 5RA > of (11.20p are confined for < f + < 1 to the 

region 



+ |9A| < ^max|i,/iocT| + (5J)2 (1.22 



Proof. See Appendix [Bj □ 
1.7.3 Asymptotic Evans function analysis 

Denoting by D{X) the "reduced" Evans function (defined Section[2|) associated with the reduced 
eigenvalue equations ()1.20p , we introduce the pair of renormalizations 

((1 - Bl/^f + 4A(^/2 + l/2(j^o)) (v+/4. + x) ^'^ 



xl/4 / n1/4 

(1 - Bl/^f + 4(^/2 + l/2a^,o)) [v+/A + 1 
1 - Bl/^JWTf + AX{^i/2v+ + l/2a^XQvl^^^'^ 



(1.23) 



and 



{{l-Bl/^JI^f + A{i,/2v+ + 

1 - Bl/^f + 4A(/x/2 + l/2a/io)^ ^'^ 
((1 - Bl/^Y + 4(^/2 + l/2a^i, 



1/4 



D{X) 



D{X) := ^D{X). (1.24) 



Intermediate behavior. 
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Theorem 1.10. On 3f?A > 0, the reduced Evans function D is analytic in A and continuous 
in all parameters except at = and Bi = y/JIov±, at which points it exhibits algebraic 
singularities (blow-up) at X = 0. The renormalized Evans functions D and D are analytic in A 
and continuous in all parameters except at {X,v^) = (0,0). 

Proof. Immediate from Propositions 12.13] and [3T6l □ 
The small-amphtude limit. 

Proposition 1.11 ([311[53]). For a > 0, > 0, > bounded, and i?J bounded away from 
y/jM)j parallel shocks are Evans stable in both full and reduced sense in the small- amplitude limit 
V-i. — > 1. Moreover, D converges uniformly on compact subsets of {3ftA > 0} as ^ 1 to a 
nonzero real constant. 

Proof. For l^;,. — = |1 — w+l sufficiently small and B^ bounded away from ^/JIq, the as- 
sociated profile must be a Lax 1- or 3-shock, whence stability follows by the small-amplitude 
results obtained by energy estimates in [31] or by asymptotic Evans function techniques in [44] . 
Convergence on compact sets follows by the argument of Proposition 4.9, |i2Sj, which likewise 
uses techniques from [44]. □ 



The large-amplitude limit. 

Theorem 1.12. For a, /xq and B^ bounded, the reduced Evans function D{X) converges uni- 
formly on compact subsets of {3iX > 0} \ {0} in the large- amplitude limit u + ^ to a limiting 
Evans function D^{X) obtained by substituting vq for v in (jl.20p . vq as in Proposition see 
Definition \3.4\ for a precise definition. Likewise, D and D converge to 

il-Bl/^o? + AX{^^/2 + l/2a^,o)f'^ 
D\X) := ^ J-^X^^D^iX) (1.25) 

(l - Bl/^)^ + 4(^/2 + l/2a/io)) 

and 

Ul-Bl/^y + 4X{^^/2 + l/2a^lo)) 

D'{X) := (^)' (1-26) 

^1 - S|/^)2 + 4(^/2 + l/2a^o) 



each continuous on 5RA > 0. Moreover, for B^ < ^/JIq, nonvanishing of on {5RA > 0} is 
necessary and nonvanishing of on {5RA > 0} is sufficient for reduced Evans stability (i.e., 
nonvanishing of D, D on {KA > 0} for Vj^ > sufficiently small. For B\ > yfJiQ, nonvanishing 
of on {3f?A > 0} is necessary and nonvanishing of on {5RA > 0} together with a certain 
sign condition on D{0) is sufficient for reduced Evans stability for ?;+ > sufficiently small^ 
This sign condition is implied in particular by nonvanishing of D(Q) on the range 

Va^<Si* <^ + max{y^,y^}. (1.27) 



D, D are real-valued for real A by construction, so that sgnl)(0) is well-defined. 
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Proof. Convergence follows by Proposition 13.61 for stability criteria, see Section 13.51 □ 

Remark 1.13. The theoretically cumbersome condition ()1.27p is in practice no restriction, 
since we check in any case the stronger condition of nonvanishing of D for 5RA > on the 

entire range y/JIo < < ^/JIq + max | \J^, \[^ } ■ 



Remark 1.14. Recall that the associated gas-dynamical shock has already been shown to 
be Evans stable for v+ > sufficiently small. Thus, not only reduced Evans stability, but full 
Evans stability, is implied for > sufficiently small by stability of the limiting function 
(resp. D^). 

Large- and small-parameter limits. 

Theorem 1.15. For a, fiQ bounded and bounded from zero, and bounded from zero, parallel 
shocks are reduced Evans stable in the limit as B\ ^ oo or ^ 0. For A bounded and 5RA > 0, 
the Evans function converges as B^ ^ oo to a constant. 

Theorem 1.16. For B* bounded, and f+ bounded from zero, parallel shocks are reduced Evans 
stable in the limit as a ^ with /io bounded, fiQ ^ with a bounded. In each case, the Evans 
function converges uniformly to zero on compact subsets of {RX > 0}, with C'^^ < \D''{X)\ < 
for Ca < |A| < C and C'^^ < I-D'^IA)] < for CfiQ < |A| < C. 

Theorem 1.17. For B^ bounded, andvj^ and iiq bounded from zero, parallel shocks are reduced 
Evans stable in the limit as a^Q — > oo. For A bounded and 5RA > 0, the Evans function D, 
appropriately renormalized, converges as crfiQ oo to the Evans function D for ()1.2ip ; more 
precisely, D ~ eCoCT^o+ci+c2A^ ^. QQ^^giant. 

Remark 1.18. Except for certain "corner points" consisting of simultaneous limits of B\ — > oo 
together with o" ^ 0, or ajjLQ oo together with — > or /Uq — > 0, our analytic results verify 
stability on all but a (large but) compact set of parameters. We conjecture that stability holds 
in these limits as well; this would be an interesting question for further investigation. 

As pointed out in ^29], the limit — > is connected with the isentropic approximation, and 
does not occur for full (nonisentropic) MHD for gas constant 7 > 1; thus, a somewhat more 
comprehensive analysis is possible in that case. Note that the reduced eigenvalue equations are 
identical in the nonisentropic case ]13^ . except with v replaced by a full (nonisentropic) gas- 
dynamical profile, from which observation the reader may check that all of the analytical results 
of this paper goes through unchanged in the nonisentropic case, since the analysis depends only 
on V, and this only through properties of monotone decrease, \vx\ < C\v\ (immediate for v 
bounded from zero), and uniform exponential convergence as x ^ ±00, that are common to 
both the isentropic and nonisentropic ideal gas cases. 

Discussion. Taken together, and along with the previous theoretical and numerical investiga- 
tions of [27] on stability of gas-dynamical shocks, our asymptotic stability results reduce the 
study of stability of parallel MHD shocks, in accordance with the general philosophy set out 
in [271 EHI [28] , rnainly (i.e., with the exception of "corner points" discussed in Remark ll.lSp 
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to investigation of the continuous and numerically well-conditioned renormalized functions D 
and D on a compact parameter-range suitable for discretization, together with investigation of 
the similarly well-conditioned limiting functions and D^. However, notice that the same 
results show that the unrenormalized Evans function D blows up as A — > 0, both in the large- 
amplitude limit f _|_ ^ and in the characteristic limits — > ^yjIov±, hence is not suitable 
for numerical testing across the entire parameter range. Indeed, in practice these singularities 
dominate behavior even rather far from the actual blow-up points, making numerical investi- 
gation infeasibly expensive if renormalization is not carried out, even for intermediate values 
of parameters/frequencies. This is a substantial difference between the current and previous 
analyses, and represents the main new difficulty that we have overcome in the present work. 

1.8 Numerical stability results 

For a given amplitude, the above analytical results truncate the computational domain to 
a compact set, thus allowing for a comprehensive numerical Evans function study patterned 
after |27l [29], which yields Evans stability in the intermediate parameter range. We then 
demonstrate Evans stability in the large-amplitude limit by (i) verifying convergence to the 
limiting Evans functions given in Theorem 11.121 (i.e., checking that convergence has occurred 
to desired tolerance at the limits of values v^, A considered), and (ii) verifying nonvanishing 
on 3f?A > of the limiting functions D^, . These computational results, together with the 
analytical results in Section 11.71 give unconditional stability for all values except for cases 
where two or more parameters blow up simultaneously as described in Remark 11.181 The 
numerical computations were performed by the authors' STABLAB package, which is written 
in MATLAB, and has been used successfully for several systems [H [?7l[29t [TTl [26l [28 ] . 

When compared to the numerical study for isentropic Navier-Stokes [U [27], this present 
system is better conditioned, yet much more computationally taxing since there are more free 
parameters to cover, i.e., (7, SJ', //q, c); the isentropic model by contrast has only two pa- 
rameters (7,w+). Since each dimension adds, roughly, an order of magnitude to the runtime, 
we upgraded our STABLAB package to allow for parallel computation via MATLAB 's paral- 
lel computing toolbox. In our main study, we computed along 30,870 semi-circular contours 
corresponding to the parameter values 

(7,u+,Si*,;Uo,cj) G [1.0,3.0] X [10~"^0.8] x [0.2,3.8] x [0.2,3.8] x [0.2,3.8]. 

In every case, the winding number was zero, thus demonstrating Evans stability; see Section [5] 
for more details. 

We also carried out a number of small studies to illustrate our analytical work in the 
limiting fixed-amplitude cases. These are briefly described below and are also given more detail 
in Section [5] 

In Figure [H we see the typical concentric structure as u+ varies on [0, 1]. Note that in the 
strong-shock limit, the output converges to the outer contour representing the Evans function 
output of the limiting system. In the small-amplitude limit, the system converges to a non-zero 
constant. Since the origin is outside of the contours, one can visually verify that the winding 
number is zero thus implying Evans stability, even in the strong-shock limit. 
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Figure 1: Renormalized Evans function output for semi-circular contour of radius 4.5 (left) 
as the amplitude varies. Parameters are = 2, fiQ = 1, a = 1, j = 5/3, with 
v+ = 10-\ 10-1^10-2, 10-2-5, 10-3, 10-3-5, 10-4, l0-'^■^10-^10-5•^ 10-6. Note the striking 
concentric structure of the contours, which converge to the outer contour in the large- amplitude 
limit (i.e., f+ 0) and to a non-zero constant in the small-amplitude limit (i.e., — > 1), in- 
dicating stability for all shock strengths since the winding numbers throughout are all zero. 
The limiting contour given by D^{X) is also displayed, but is essentially identical to nearby 
contours. When the image is zoomed in near the origin (right), which is marked by a crosshair, 
we see that the curves are well behaved and distinct from the origin. Also clearly visible is the 
theoretically predicted square-root singularity at the origin of the limiting contour, as indicated 
by a right angle in the curve at the image of the origin on the real axis. 

In Figure El we illustrate the convergence of the Evans function as -B* oo. Note that the 
contours converge to zero, but they are stable for all finite values of -B*. Stability is proven 
analytically in Theorem 11.151 by a tracking argument. Prior to this computation, however, a 
significant effort was made to prove stability with energy estimates, but these efforts were in 
vain since the Evans function converges to zero as oo. 

In Figure El we see the structure as 0. Once normalized (right), we see that the 

structure is essentially unchanged despite a large variation in /ig; in particular, the shock layers 
are stable in the /xq — > limit. This was proven analytically in Theorem 11.161 

Finally, in Figure HI we see the behavior of the Evans function in the case that r = ^/{2^ + 
r/) ^ oo. This is the opposite case of that considered in |l3j. As we show in Proposition 14.31 this 
case can be computed by disengaging the shooting algorithm and just taking the determinant 
of initializing e-bases at ±00. Notice that in this limit the shock layers are also stable. 

1.9 Discussion and open problems 

Our numerical and analytical investigations suggest strongly (and in some cases rigorously 
prove) reduced Evans stability of parallel ideal isentropic MHD shock layers, independent of am- 
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Figure 2: Evans function output for semi-circular contour of radius 5 (left) and a zoom-in of 
the same image near the origin (right). Parameters are V-^- = 10~^, = 1, a = 1, ^ = 5/3, 
with Bl = 2,3.5,5,10,15,20,25,30,35,40. Note that the contours converge to zero, which is 
marked by a cross hair, as i?* ^ oo. 

plitude, viscosity and other transport parameters, or magnetic field, for gas constant 7 G [1,3], 
indicating that they are stable whenever the associated gas-dynamical shock layer is stable. 
Together with previous investigations of [27| indicating unconditional stability of isentropic 
gas-dynamical shock layers for 7 S [1,3], this suggests unconditional stability of parallel isen- 
tropic MHD shocks for gas constant 7 G [1,3], the first such comprehensive result for shock 
layers in MHD. 

It is remarkable that, despite the complexity of solution structure and shock types occurring 
as magnetic field and other parameters vary, we are able to carry out a uniform numerical Evans 
function analysis across almost (see Remark ll.lSp the entire parameter range: a testimony 
to the power of the Evans function formulation. Interesting aspects of the present analysis 
beyond what has been done in the study of gas dynamical shocks in \27\ [29] are the presence of 
branch singularities on certain parameter boundaries, necessitating renormalization of the Evans 
function to remove blow-up singularities, and the essential use of winding number computations 
on Riemann surfaces in order to establish stability in the large-amplitude limit. The latter 
possibility was suggested in [16] (see Remark 3, Section 2.1), but to our knowledge has not up 
to now been carried out. 

We note that Freistiihler and Trakhinin jj3] have previously established spectral stability 
of parallel viscous MHD shocks using energy estimates in the regime 

r := fi/ilfi + T]) < 1, 

whenever < 2yJ^QV^ (translating their results to our setting s = —1), which includes all 
1- and intermediate-shocks, and some slow shocks {B\ > yJJiQvZ). Recall that we have here 
followed the standard physical prescription r/ = — 2/x/3, so that /x/(2/x + ff) = 3/4, outside the 
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Figure 3: Evans function output for semi-circular contour of radius 5 (left) together with a 
renormalized version of the contours (right), where each contour is divided by its rightmost 
value, thus putting all contours through z = 1 on the right side. Although these results 
are typical, the parameters in this example are Bl = 2, f+ = 10~^, cr = 1, 7 = 5/3, with 
= 10--5, 10-\ 10-1-5, 10-2, 10-2■^ 10-3, 10-3-5, 10-4, 10-4-5, 10-5. Note that the renormal- 
ized contours are nearly identical. This provides a striking indication of stability for all values 
of ^0 in our range of consideration, and in particular for fiQ 0. 

regime studied in [13j. Thus, the two analyses are complementary. It would be an interesting 
mathematical question to investigate stability for general ratios fi/{2fi + ij). See Section [4.41 for 
further discussion of this issue. Here we study only the limit r ^ 00 complementary to that 
studied by [13], the case r = 3/4 suggested by nonmagnetic gas dynamics, and the remaining 
cases in the r — > limit left open in [13]. Other r- values may be studied numerically, but were 
not checked. 

Stability of general (not necessarily parallel) fast shocks in the small magnetic field limit 
has been established in [19] by convergence of the Evans function to the gas-dynamical limit, 
assuming that the limiting gas-dynamical shock is stable, as has been numerically verified for 
ideal gas dynamics in [271 [29l I28j . Stability of more general, non-gas-dynamical shocks with 
large magnetic field, is a very interesting open question. In particular, as noted in [|8j, one- 
dimensional instability, by stability index considerations, would for an ideal gas equation of 
state imply the interesting phenomenon of Hopf bifurcation to time-periodic, or "galloping" 
behavior at the transition to instability. For analyses of the related inviscid stability problem, 
see, e.g., dHUHHO] and references therein. 

Another interesting direction for further investigation would be a corresponding compre- 
hensive study of multi-dimensional stability of parallel MHD shock layers, as carried out for 
gas-dynamical shocks in [28j. As pointed out in ^13j, instability results of [Gj |49] for the corre- 
sponding inviscid problem imply that parallel shock layers become multi-dimensionally unstable 
for large enough magnetic field, by the general result [56^ [5T] that inviscid stability is necessary 
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Figure 4: Renormalized Evans function in the r = oo case. Parameters are f+ = 

10^^, 10^^, 10^^, 10^^, /io = 1, o" = 1, 7 = 5/3. We also have a semi-circular radius of 4.5 with 

B\ = 2 (left), and a semi-circular radius of 1 with = 0.5 (right). 

for viscous stability, so that in multi-dimensions instability definitely occurs. The question in 
this case is whether viscous effects can hasten the onset of instability, that is, whether viscous 
instability can occur in the presence of inviscid stability. 

2 The Evans function and its properties 

We begin by constructing carefully the Evans function associated with reduced system ()1.20p . 
and recalling its basic properties for our later analysis. 

2.1 The Evans system 

The reduced eigenvalue equations (11.201) may be written as a first-order system 
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(2.1) 



or 

W' = A{x,X)W, (2.2) 

indexed by the three parameters B^, a, and /Iq := afiQ. Recall that we have already fixed /i = 1 
and (2/i + tj) = 4/3. 
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2.2 Limiting subspaces 

Denote by 



/ 1/^ \ 

\v± v±/iJ, —aB\v± 
criJ,ov± 

\ -Blv±/fi Xv± cT//of± J 

the limiting coefficient matrices associated with (I2.ip - ()2.2p . 



A±{X) := hm ^(x,A) 

a;— >±cx3 



(2.3) 



Lemma 2.1. For 5RA > 0, A 7^ 0, and 1 > > 0, each of A± has two eigenvalues with 
strictly positive real part and two eigenvalues with strictly negative real part, hence their stable 
and unstable subspaces S± and U± vary smoothly in all parameters and analytically in A. For 
w+ > they extend continuously to A = 0, and analytically everywhere except at B\ = yJ^ov±, 

where they depend smoothly on {1 — B^/ iiqv±Y + 4Ac, c := ^ + ^^^^^i ^ • 

Proof. By standard hyperbohc-parabolic theory (e.g., Lemma 2.21, [51]), A± have no pure 
imaginary eigenvalues for 3f?A > 0, A 7^ 0, for any Vj^ > Q and parameter values a, /io, -B^, 
whence the numbers of stable (negative real part) and unstable (positive real part) eigenvalues 
are constant on this set. By homotopy taking A to positive real infinity, we find readily that 
there must be two of each. 

Alternatively, we may see this directly by looking at the corresponding second-order sym- 
metrizable hyperbolic-parabolic system 

w\ , / 1 —Bi/^qv±\ f w\ _ f fJ-/v± \ fw 



X + 



a) \—B\ 1 y \'^/ X V l/o'Aio^^±y \a 

and applying the standard theory here, specifically, noting that eigenvalues consist of solutions 

^-(-"U: (2.4, 

so that /i = i/c by a straightforward energy estimate yields JRA < —Ok"^ for 6 > Q. 

Applying to reduced system (j2.4p Lemma 6.1 [38] or Proposition 2.1, [55j . we find further 
that, whenever the convection matrices 

&:=(_^. -''''r') (2.5) 

are noncharacteristic in the sense that their eigenvalues a = 1 it , ^ are nonzero, these 

° Vw± ' 

subspaces extend analytically to A = 0. 

Finally, we consider the degenerate case that B^ = y/JiQV± and the convection matrix fi is 
characteristic. Considering (|2.4p as determining A//U as a function of /i for /i small, we obtain, 
diagonalizing B and applying standard matrix perturbation theory ^34] that in the nonzero 
eigendirection rj of /3, associated with eigenvalue /3j / 0, A//i ~ and, inverting, we find that 
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/ij(A) extends analytically to A = 0. In the zero eigendirection, on the other hand, associated 
with left and right eigenvectors / = (1/2, 1/2B\)^ and r = {B, I)-'", we find that 

A ~ -//(I - BI/^hqv±) + c^^ + . . . , 

where 

leading after inversion to the claimed square-root singularity. Likewise, the stable eigendirec- 
tions of A± associated with vary continuously with A, converging for B\ = ^JJIqv± and A = 

to {w, a, w' ,a')^ = , where r is the zero eigendirection of j3. This accounts for three eigen- 
values /i lying near zero, bifurcating from the three-dimensional kernel of A± near a degenerate, 
characteristic, value of The fourth eigenvalue is far from zero and so varies analytically in 
all parameters about A = 0. □ 

Remark 2.2. Remarkably, even though the shock changes type upon passage through the points 
B\ = ^/JM)V±, the stable and unstable subspaces of A± vary continuously, with stable and 
unstable eigendirections coalescing in the characteristic mode. 



2.3 Limiting eigenbases and Kato's ODE 

Denote by n_|_ and n_ the eigenprojections of yl+ onto its stable subspace and onto its 
unstable subspace, with A± defined as in (j2.3p . By Lemma 12.11 these are analytic in A for 
5RA > 0, v+ > 0, except for square-root singularities at A = for B\ = ^/JIov±. Introduce the 
complex ODE [S] 

i?' = n'i?, R{\q) = Ro, (2.6) 

where ' denotes d/dX, Aq is fixed with 3?Ao > 0, 11 = II-i-, and i? is a 4 x 2 complex matrix. By a 
partition of unity argument [34| , there exists a choice of initializing matrices Rq that is smooth 
in the suppressed parameters v+, a, and fiQ, is full rank, and satisfies n(Ao)i?o = -^o; that 
is, its columns are a basis for the stable (resp. unstable) subspace of A^ (resp. A_). 

Lemma 2.3 ([Ml [53j). There exists a global solution R of (j2.6p on {5RA > 0}, analytic in A 
and smooth in parameters V-^- > 0, B^ > 0, a > 0, and fiQ > except at the singular values 
A = 0, B^ = ^HoV±, such that (i) rankii = ranki?*^, (ii) HR = R, and (Hi) liR! = 0. 

Proof. As a linear ODE with analytic coefficients, (j2.6p possesses an analytic solution in a 
neighborhood of Aq, that may be extended globally along any curve, whence, by the principle of 
analytic continuation, it possesses a global analytic solution on any simply connected domain 
containing Aq [34] • Property (i) follows likewise by the fact that R satisfies a linear ODE. 
Differentiating the identity 11^ = 11 following [M] yields 1111' -|- II'II = 11', whence, multiplying 
on the right by 11, we find the key property 

nn'n = o. (2.7) 
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From (j2.7p . we obtain 

{HR - R)' = {H'R + HR! - R!) = H'R + (H - 1)11' R = HH'i?, 
which, by nil'n = and 11^ = 11 gives 

(Hi? - R)' = -nn'(ni? - r), {hr - r){\o) = o, 

from which (ii) fohows by uniqueness of solutions of hnear ODE. Expanding liR' = R and 
using IiR = R and nn'H = 0, we obtain IIR' = lili'IiR = 0, verifying (iii). □ 

Remark 2.4. Property (iii) indicates that the Kato basis is an optimal choice in the sense that 
it involves minimal variation in R. It is useful also for computing the Kato basis in different 
ways f3^. see Appendix O 



2.4 Characteristic values: the regularized Kato basis 

We next investigate the behavior of the Kato basis near A = and the degenerate points 
= y/jM)V± at which the reduced convection matrix I3± of (|2.5p becomes characteristic in a 
single eigendirection. 

Example 2.5. A model for this situation is the eigenvalue equation for a scalar convected heat 
equation Xu+rju' = u" with convection coefficient rj passing through zero. The coefficient matrix 
for the associated first-order system is 

A:=f» 'V ,2.8) 



As computed in Appendixl^ the stable eigenvector of A determined by Kato's ODE (j2.6p is 

(r/V4 + 1)1/4 



Riv, A) := l^^y^^^ji/. (1' - VWITA j , (2.9) 

which, apart from the divergent factor ^^jyj^^jT/l > smooth function of Y^r/2/4 + A. 

The computation of Example (|2.5p indicates that the Kato basis blows up at A = as ((1 — 

Bl / ^/JM)V±)'^ +4\) as B'l crosses characteristic points ^/iof± across which the shock changes 
type, hence does not give a choice that is continuous across the entire range of shock profiles. 
However, the same example shows that there is a different choice (1, —ri/2 — ^ rf' /A + A )'^ that 
is continuous, possessing only a square-root singularity. We can effectively exchange one for 
another, by rescaling the Kato basis as we now describe. 

Following [HITH], associate with bases R± = (i?i,i?2)^ tlie wedge (i.e., exterior algebraic) 
product TZ± := {Ri A i?2)=^. 
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Lemma 2.6. The "regularized Kato products" 

1/4 



((1 - Bl/^/JWrf + AX{^l/2v+ + l/2a^,ov] 
((1 - 5J/^w^)2 + 4(/x/2r;+ + l/2a^lQ^ 



7^+ := rT7r« ^ ^2^) (2.10) 



and 



1/4 



((1 - 5iV^)2 + 4A(^/2 + l/2a^o 

^- ■■=- '^{Rt^Rt) (2-11) 

:i - Bl/^f + 4W2 + l/2a/io) 



are analytic in A and smooth in remaining parameters on all of X > 0, V-^- > 0, B'l > 0, a > 0, 
/io > except the points A = 0, i?* = ^/iof±, where they are continuous with a square-root 

singularity, depending smoothly on ^{1 — B^/ ^/^o^^'i)^ + 4A. Moreover, they are bounded from 
zero (full rank) on the entire parameter range. 

Proof. A computation like that of Example 12.51 applied to system (j2.3p , replacing r/ with the 
characteristic speed 1 — B^/ ^ ^qv± and introducing a diffusion coefficient c± = \ + ^^j^^ij ^ , 
i.e., considering 

\u + riu' = cu", (2.12) 

shows that, for an appropriate choice of initializing basis Rq in (j2.6|) . there is blowup as 
{\,B\) —> (0, y^^lol^) at rate 



((l-Si7Vw^)' + 4Ac±) 



-1/4 



in a basis vector involving the characteristic mode, while the second basis vector remains 
bounded and analytic, whence the result follows. For the derivation of approximate equation 
(j2.12p . see the proof of Lemma |2.1[ □ 

Remark 2.7. A review of the argument shows that estimate (j2.10p derived for fixed vj^ remains 
valid so long as |A| < C\^\ <^ v\. Different asymptotics hold for vj^ < Cy^|A[; see Section [Ql 

Remark 2.8. Lemma \2.6\ (by uniform full rank) includes the information that the unregularized 
Kato bases blow up at rate A^"*^/^ at B^ = ^/^Iov±. 

2.5 Conjugation to constant-coefficients 

We now recall the conjugation lemma of [41j. Consider a general first-order system 

W = A{x,X,p)W (2.13) 

with asymptotic limits A± as x ^ ±oo, where p £ denotes up-to-now-supressed model 
parameters. 
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Lemma 2.9 (dUHl]). Suppose for fixed 9 > and C > that 

1^- A±|(x,A,p) < C7e-^l^l (2.14) 

for X ^ uniformly for {X,p) in a neighborhood of {Xo,po) and that A varies analytically 
in A and smoothly (resp. continuously) in p as a function into L°^{x). Then, there exist in 
a neighborhood of (Ao,Po) invertible linear transformations P+{x, X,p) = I + 0+(x,A,p) and 
P-{x, X,p) = /+0_(x, X,p) defined on x >0 and x < 0, respectively, analytic in X and smooth 
(resp. continuous) in p as functions into L°°[0, ±00), such that 

|e±| < C7ie-^>l for x^O, (2.15) 
for any < 9 < 9, some Ci = Ci{9,9) > 0, and the change of coordinates W =: P±Z reduces 

(iron to 

Z' = A±Z for X ^ 0. (2.16) 

Proof. The conjugators P± are constructed by a fixed point argument [41j as the solution of an 
integral equation corresponding to the homological equation 

P' = AP-A±P. (2.17) 

The exponential decay (j2.14p is needed to make the integral equation contractive in L°°[M, +00) 
for M sufficiently large. Continuity of P± with respect to p (resp. analyticity with respect to A) 
then follow by continuous (resp. analytic) dependence on parameters of fixed point solutions. 
Here, we are using also the fact that (j2.14p plus continuity of A from p L°° together imply 
continuity of e^l^l(j4 — A±) from p into L°^[0,ibcxD) for any < 9 < 6, in order to obtain the 
needed continuity from p L°° of the fixed point mapping. See also [lll[20]. □ 

Remark 2.10. In the special case that A is block- diagonal or -triangular, the conjugators P± 
may evidently be taken block- diagonal or triangular as well, by carrying out the same fixed-point 
argument on the invariant subspace of ()2.17p consisting of matrices with this special form. This 
can be of use in problems with multiple scales; see, for example, the proof in Section^of Theorem 

irm (a ^ 0;. 

2.6 Construction of the Evans function 

Definition 2.11 ( [Ml ED [52 J ) . The Evans function is defined on I^M) on 5RA > 0, v+ > 0, 

> 0, a > 0, Ho > as 

DiX,p) := det{P+Rt,P+Rt,P-Ri,P-R2)\.=o 
= {P+Rt AP^R^ AP-R^ AP-R2\x=o), 

where {■) of a full wedge product denotes its coordinatization in the standard (single- element) 
basis ei A 62 A 63 A 64, where Cj are the standard Euclidean basis elements in C^. 
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Definition 2.12. The regularized Evans function is defined as 

((1 - i3iV^)2 + 4A (/x/2 + l/2a//o)) 
D{X,p) := ^ 



((1 - + 4A(/i/2i;+ + 1/2(7^0^^^ 



1/4 (2.19) 
D{X,p) 



((1 - i3f/V7w^)2 + 4{fi/2v+ + l/2afiovl)) ^'^ 

= AP_i?„U=o), 

where V±{Ri A R2) ■= P±Ri f\P±R2 denotes the "lifting" to wedge product space of conjugator 
P±- 

Proposition 2.13. The Evans function D is analytic in A and smooth in remaining parameters 
on all of X > 0, v+ > 0, Bl > 0, a > 0, Ho > except the points X = 0, Bf = ^Hov±, where it 
blows up as 

\l - Bl/,/JI^f + AX{ii/2v± + l/2c7^o4)^ 



The regularized Evans function D is analytic in X and smooth in remaining parameters on the 
same domain, and continuous with a square-root singularity at (A, B\) = (0, ^/Jiov±), depending 



smoothly on \J — Bl/ ^ hqv±Y + 4A. 



Proof. Local existence/regularity is immediate, by Lemmas 12.31 12.6] and 12.91 Proposition 11.1^ 
and Definitions [2lT| [2J2l Global existence/regularity then follow [38, 441 EH [52] by the 
observation that the Evans function is independent of the choice of conjugators P± (in general 
nonunique) on the region where A± are hyperbolic (have no center subspace), in this case 
{^X > 0} \ {0}E □ 

Remark 2.14. Evidently, for B^ 7^ \/fJ'OV±, Evans stability, defined as nonvanishing of D on 
5RA > is equivalent to nonvanishing of the regularized Evans function D on ?R.X > 0. On 
the other hand, D is continuous throughout the physical parameter range, making possible a 
numerical verification of nonvanishing, even up to the characteristic points B^ = y/JIov±. 

Remark 2.15. An alternative, simpler and more general regularization of the Evans function 
is 

where and \V-'R.-\ denote norms of {VTZ)± in the standard basis CiAcj, i / j. Though 

not analytic, it is still C°° wherever D is analytic, and its zeros agree in location and multiplic- 
ity with those of D, and is continuous wherever the stable (unstable) subspaces of (^-) f^i^^ 



In Evans function terminology, the "region of consistent splitting" [T] 1161 1511 152) . 
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continuous. Indeed, it is somewhat more faithful than the usual Evans function to the original 
idea fj^ of a quantity measuring the angle between subspaces. The disadvantage of this regular- 
ization is that it eliminates structure (analyticity, asymptotic behavior) that has proved quite 
useful both in verifying code by benchmarks, and in interpreting behavior /trends l27[\ll\\29^ \28 ^ . 

Proposition 2.16. On {SRA > 0} \ {0}, the zeros of D (resp. D) agree in location and 
multiplicity with eigenvalues of L. 

Proof. Agreement in location- the part tliat concerns us here- is an immediate consequence 
of the construction. Agreement in multiplicity was established in |14^ I15j. For an alternative 
argument, see [55l|38]. □ 

3 The strong shock hmit 

We now investigate behavior of the Evans function in the strong shock limit f + — > 0. By Lemma 
12.91 Proposition ll.il and Corollary ll.2l this reduces to the problem of finding the limiting Kato 
basis i?+ at +oo as f+ 0. That is, this is a "regular perturbation" problem in the sense 
of [44 1 l29], and not a singular perturbation as in the much more difficult treatment of the 
gas-dynamical part {v,ui) done in [27] . On the other hand, we face new difficulties associated 
with vanishing of the limiting Evans function at A = and branch points in both limiting and 
finite Kato flows, which require additional stability index and Riemann surface computations 
to complete the analysis. 

3.1 Limiting eigenbasis at +oo as f+ ^ 0, |A| > 6* > 

Fixing n = I without loss of generality, we examine the limit of the stable subspace as ?;_)_—> 
of 



/ 


1 









v+ 





-aBlv+ 











cr^o 


V 




Xvl 





(3.1) 



Making the "balancing" transformation 

A+ := v-^/^TA+T-\ T:=diag{t;i/', 1,1,1} (3.2) 
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and expanding in powers of ^4. , we obtain 



+ V 



/O 1 0\ 

A 



Vo 0/ 

/O 

3/2 



, 1/2 

+ f+ 








\ 

1 -aB 
cj/io 
-Sf A / 



'1 



(3.3) 



\ 




\0 cj/io/ 



Noting that the upper lefthand 2x2 block of has eigenvalues ib\/A bounded from zero, 
we find [M] that ^+ has invariant projections Hi = 7^i£| and 112 = '^2^2 within 0(f_j_ ) 
of the standard Euclidean projections onto the first-second and the third-fourth coordinate 
directions, i.e., 

/I 0\ /o 0\ 



7^l 



1 

VO 0/ 



+ 0(7;^'), 7^2 





1 
VO 1/ 



o« ), 



r / 1 0^ ^n^ l/2^ ^ /O 1 

^i = 'o 1 oj+^("+ )'^2=lo 

1/2 



Indeed, looking more closely- expanding in powers of vj_ and matching terms- we find 
after a brief calculation 



7^2 



Looking at L\Aj^TI, 
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D* 1/2 \ 
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vo 


1 / 






7^l 







+ 0{v+), £2 



,1/2 



Blv 



1 0' 
1, 



+ 0{v+). 



1/2 



+ 0('u+) and noting that ibV^ are spectrally separated 



by the assumption |A| > > 0, we find that the stable eigenvector within this space is 

(1, - Vt'+/4 + Af + 0{vf), 
and thus the corresponding stable eigenvector within the full space is 



i?i = (1, - ^v+/A + A, 0, 0)^ + 0( 



,l/2x 
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— 1/2 

Looking at C2^+T^2 



'7^0 \ 3/2n noting that the eigenvalues —aiiQVA./2± 

a'^n'^v'^/A + anoX of the principal part are again spectrally separated so long as afiQ > are 
held fixed, we find that the stable eigenvector within this space is 

T 



{l,-v+/2-^vl/A + X/aiso) +0(t;f), 
and thus the corresponding stable eigenvector within the full space is 



R2 



O{v]!^),0,l,-v+/2 



vl/A + X/afio ) +0{v^f). 



Converting back to original coordinates, we find stable eigendirections T ^Ri = (1, 0, 0, 0)"^ + 



0{v_ 



1/2^ 



and 



or, using an appropriate linear combination 








*, 0, 1, -V+/2 - + X/afio ) + 0{v^^), 



( 



+ 0(1-, 



1/2^ 



Ri 






1 



\-v^l2 - ^vl/A + X/aixo) 



+ 0{v. 



1/2. 



(3.4) 



Finally, we deduce the limiting Kato ODE flow as u+ — > 0. A straightforward property of 
the Kato ODE is that it is invariant under constant coordinate transformations such as (j3.2p . 
Thus, we find, for appropriate initialization, that i?5 = (r, 0,0, 0)^, where (r, s)-^ is the Kato 



eigenvector associated with 





A V 



1 



/2 , or (by the calculation of Example 12. 5^ setting rj = v. 



1/2, 



/f+/4+ 

U+/4 + A 



:i, 0,0,0)^ 



(3.5) 



with limit 

i?o = (A-l/^ 0,0,0)^. (3.6) 

Similar considerations yield a second limiting solution R2 = (0, 0, r, s)^, where (r, s)"^ is the 

'0 afiQ 



Kato eigenvector associated with 



A afiov+ 



or 



R2 



vl/i + X/ano^ 



(0, 0, 1, -V+/2 - ^vi/A + X/^o? 



with limit 



i?o=(o,0,A-V4,-AVVV^)) . 
We collect these observations as the following lemma. 



(3.7) 
(3.8) 
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Lemma 3.1. On compact subsets of {JRA > 0} \ {0}, Ri and R2 converge uniformly in relative 
error to fixed (i.e., independent of X) linear combinations of Ri and R2 as defined in (j3.6p and 

Remark 3.2. The above computations show that the formulae for R^ remain valid so long as 
|A| ^ v'^. Recall, for |A| <^ v'^, the behavior is as described in Lemma \2. 61 This leaves only the 
case \X\ v\ unexamined. 

3.2 Limiting behavior at +00 as 11+, A ^ 

As suggested by the different behavior for |A| ^ f ^ and |A| ^ i^^, behavior in the transition 
zone |A| ~ appears to be rather comphcated, and so we do not attempt to describe either 
the hmiting subspace or hmiting Kato flow as w+ and A simultaneously go to zero, recording 
only the following topological information. 

Lemma 3.3. In rescaled coordinates {w,w' ,a,a' /v), for f + > sufficiently small, the Kato 
product R^ A i?^ defined above is analytic for 5RA > —6, 9 > sufficiently small, except at two 
(possibly coinciding) singularities Xi, X2 near the origin, each of fourth-root type and blowing 
up as (A — Aj)~^/^. 

Proof. Equivalently, by the computation of Example 12.51 we must show that each of the stable 
eigenvalues ai, 02 of collide with unstable eigenvalues at precisely one point Aj, which is a 
branch point of degree two. Computing the characteristic polynomial p{X, a) := det(yl+ (A) — a) 
with the aid of (|2.4p . we obtain j)( A. a) = {a'^—v+a—Xv+){a'^—aiJ.ov'^a—Xafj,ov'^)—a{Bl)'^v'^a'^, 
a quadratic in A. Taking the resultant of p with d^p, we therefore obtain a quadratic polynomial 
q{X) whose roots Aj are the points at which ^+(A) has double eigenvalues. Noting that Ai = 
A2 = for = 0, we find by continuity that they lie near the origin for sufficiently small. 

Noting that d^p = 8a — 2{v+ + aiiov'^), we find that Ai = A2 only if q = (l/4)(z;+ + afj,ov'^). 
Plugging this into the linear equation d'^p{X, a) = in A gives the further information (2t;+ + 
0{vl))X = u^(-l/4 - a{Blf) + 0{vl), hence A ~ v+{-l/8 - a{Blf /2) > v\ for v+ small. 
But, in this case, the analysis of 13.11 implies that this coalescence represents a pair of branch 
points of degree two and not a single branch point of degree four; see Remark 13.21 The same 
analysis prohibits the possibility that either of Xj represents a branch point of degree four, 
hence they must each be degree two or three. Finally, the global behavior described in Lemma 
13.11 excludes the possibility that they be degree three, leaving the asserted result as the only 
possible outcome. □ 



3.3 Limiting subspaces at —00 as A ^ 

Case {i){\B\\/^ > 1) For fi =l,v- = l, 1^ becomes 
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whence we find by a standard limiting analysis [S5l [38] as A — > that the unstable subspace 
of A^, expressed in coordinates (to, a, u)', a'), is spanned by the direct sum of = 



and 



Rr, = i 1 , where ri is the stable subspace of 
subspace of 

1 
^0 1/cr/io 

with fi2 the associated eigenvalue. 
By direct computation, ri = (1 



-Bl 



1 



-Bl 



1 



1 -Bl/fio 



and S2 is the unstable 
(3.10) 



fio)^ , while 



S2 



and 



^ (1 - g/uo) - ^(1 - a/io)2 + 4a(;& 

-2(7/io5t 



*)2x T 



1, 



M2 



(1 - afio) + y/il-afioy + 4a{Biy 
1 + afio) + y/{l-afio? + Aa{Biy 



from which we recover expressions in standard coordinates [w^w' , a, ) of 



and 



Rr, 



1,0, -VW),0 



1 



V" 



(l+CTAto)+V(l-aw)^+4a(g*)2 

2 

(l-aAto)-v/(l-'^w)^+4t7(gn^ 

p^N (l+<TMo)+V(l-aw)"+4a(i?*)- 



(3.11) 



(3.12) 



Case (ii)(|i?J|/iy/io < 1) In this case, the unstable subspace of is spanned by the direct 
sum of (si,/iisi)-^ and {s2,^J^2S2), where Sj, fij are the unstable eigenvectors, eigenvalues of 
()3.10p . giving, by a similar computation as above, 



R-, 



2 



(1-(7Mo)+V(1-^/^o)^+4(t(B?? 



-{2afioB 



(i+<7mo)-V'(i-'^w)^+4^CbF? 

1 ) {l-af,o)-y/(l-af,or+MB-*y J 



(3.13) 



and i?2 as in (j3.12p . 
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Remark 3.4. The precise form of the eigenbases is not really important here, only the fact that 
in case (i) there is a limiting direction (j3.1ip corresponding to a nondecaying, zero- eigenvalue 
mode, whereas in case (ii) all solutions asymptotic to Span{i?^, i?^} decay exponentially as 
X — > — oo. 

3.4 The limiting Evans function 

Definition 3.5. We define the limiting Evans function as the Evans function associated with 
the limiting ODE (j2.ip with v = ■O'', "0*^ as defined in Corollarv \l.^ with (indeterminate for 
this system, since is almost empty) taken as := lim^^_^o-R+ computed above in ()3.1ip . 
(|3.12p . and the renormalizations D^, as in ()1.25p . ()1.26p . 

Proposition 3.6. Appropriately normalized^ D D^, D , and D uniformly on 

compact subsets o/{5RA > 0} \ {0}, up to a constant factor independent of X. Moreover, is 
continuous on > 0} and analytic except for a square-root singularity at X = 0. Both D and 
extend meromorphically to B{0,r), forr, v+ > sufficiently small, with a single square- 
root singularity A~^/^ at the origin, and D with a pair of fourth-root singularities (A — Xi)~^^^ 
and (A — X2)~^^^ for Xj G B{0,r), with D on dB{0,r) for these extensions as well. 

Proof. Convergence of D on {3f?A > 0}\{0} follows by Lemmas 12.91 and 13 . 1 \ Proposition ll.il and 
Corollary 11.21 whereupon convergence of D and D follows by comparison of (jl.23p and (jl.25p 
and of (jl.24p and (11.261) . Regularity of follows by Lemma 12.91 and regularity of formulae 
()3.6p . (j3.8p . as does holomorphic extension to B{0,r). Holomorphic extension of D and the 
asserted description of singularities follows by Lemmas 12.91 and 13.31 □ 

3.4.1 Behavior near A = 

At the origin, we have the following striking bifurcation in behavior of D^. 
Lemma 3.7. For Bl > D^{0) = 0. For < Bl < D°{0) / 0. 

Proof. The first assertion follows from the fact that, by (13. lip , for B^ > ^JJIq, both the initializ- 
ing eigenvector = (1, 0, —^/JIq, 0)^ of and the initializing eigenvectors R^, i?2 +oo are 
preserved by the flow of ()2.ip when A = 0, for any value of v^, corresponding to the fact that 
constant w = wq, a = oq are always solutions of (jl.20p when A = 0. Thus, for Bi > -^^uq, the 
first, third, and fourth columns in the determinant (|2.18p defining D^, consist of multiples of 
(1, 0, — y7*o, 0)"^, (1, 0, 0, 0), and (0, 0, 1, 0), hence the determinant is zero. The second assertion 
follows similarly from the observation that for Bl < ^/JIq, the solutions of (j2.ip corresponding 
to i?^, R2 at A = are exponentially decaying as x — > —00, hence independent of the constant 
solutions corresponding to the initializing eigenvectors Ri, R2 at +00. □ 

Remark 3.8. As the proof indicates, the bifurcation described in Lemma \3. 7| originates in the 
nature (i.e., decaying vs. constant) of solutions as x ^ —00, corresponding to change in type 

^As done automatically by our method of numerical initialization; see Section [5] 
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of the underlying shock. Generically we expect that vanishes to square-root order at X = 
for B\ > ^/JIq, since it has a square-root singularity there. 

3.5 Proof of the limiting stability criteria 

Proof of Theorem \1.1SX By Theorem 11.91 Proposition 13.61 s-^icl properties of limits of analytic 
functions, it suffices to consider the case that |A| and |f+| are arbitrarily small. Denote the 
Evans function for a given as 0""+ , suppressing other parameters. By Lemmas 12.91 and 13.31 
we may for u+ sufficiently small extend meromorphically to a ball i3(0,r) about A = 0, 
and the resulting extension is analytic (multi- valued) except at a pair of branch singularities Ai 
and A2 at which behaves as dj{X — Aj)^^/^ for complex constants dj. Making a branch cut 
on the segment between Ai and A2 as in Figure [5l we may view 0""+ as an analytic function on 
a slit, two-sheeted Riemann surface obtained by circling the deleted segment A1A2. Applying 
Proposition 13.61 again, we find that 

0"^+ {X) ^ D^{X) ^ ca\-^'^ + ci 

on dB{0,r) as — > 0, where Cj are complex constants. 

By Lemma 13. 7| cq 7^ for Bl < y/JIo. For B\ > y^^o, cq = 0, and the condition that 

~ not vanish at the origin is the condition that ci 7^ 0. Taking the winding number 
of around dB{0,r), therefore, on the two-sheeted Riemann surface we have constructed- 
that is, circling twice as D'"+ varies meromorphically- we obtain in the first place winding 
number negative one, and in the second (assuming ci 7^ 0) winding number zero. Subtracting 
the winding number about the segment A1A2, necessarily greater than or equal to negative one 
by the asymptotics of D'"+ at Aj, we find by Cauchy's Theorem/Principle of the Argument that 
for B^ < there are no zeros of within B{0, r) \ A1A2, concluding the proof in this case. 

For Bl > ^/JIo, we find that there is at most one zero of within B{0,r) \ AiA2. To 
complete the proof, we appeal as in [11] to the mod- two stability index of [HI [501 EI] > which 
counts the parity of the number of unstable eigenvalues according to its sign, and is given by 
a nonzero real multiple of D'"+{0). To establish the theorem, it suffices to prove then that this 
stability index does not change sign, since we could then conclude stability by homotopy to a 
limiting stable case a— >+ooori?J — >+cxd. (Alternatively, we could check the sign by explicit 
computation, but we do not need to do so.) Recall that D'"+{0) is a nonvanishing real multiple 
of the product of the hyperbolic stability determinant and a transversality coefficient vanishing 
if and only if the traveling wave connection is not transverse. 

As noted already in Proposition 11.41 the hyperbolic stability determinant does not vanish 
for Lax 3-shocks, so is nonvanishing for Bl > y/JiQ. The transversality coefficient is an Evans 
function-like Wronskian of decaying solutions of the linearized traveling-wave ODE 

hence converges by Lemma 12.91 to the corresponding Wronskian for the limiting system with 
■i) replaced by . But, this limit must be nonzero wherever ci is nonzero, or else would 



29 




Figure 5: Winding number on two-sheeted Riemann surface. 



vanish at A = to at least order A due to a second hnear dependence in decaying as well as 
asymptotically constant modes, and so ci = in contradiction to our assumptions. Therefore, 

transversality holds by assumption for < i?* — ^/JIq < max{ ^J^, } ^+ sufficiently 
small. 

On the other hand, an energy estimate like that of Section [B . 31 sharpened by the observation 
that \vl\ < V, improving the general estimate Vx < 7^, yields transversality of ()3.14p for 

-^1 ~ ^ maxjy'^, }• Thus, we have transversality for all B\ > ^/JIq, and we may 
conclude by homotopy to the stable -B* — > oo limit that the transversality coefficient has a sign 
consistent with stability, that is, there are an even number of nonstable zeros ?R.X > of the 
Evans function D'"+ for V-^- sufficiently small. Since we have already established that there is 
at most one nonstable zeros of D'"+ , this implies that there are no nonstable zeros, yielding 
stability as claimed. □ 



Remark 3.9. From Lemma \3. 7\ there might appear to be inherent numerical difficulty in ver- 
ifying nonvanishing of D for B\ less than but close to ^/JIq, since D vanishes at the origin for 
^1 — \/A*o- However, this is only apparent, since we know analytically that does not vanish 
at, hence also near, A = for B\ < 
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4 Further asymptotic limits 



In this section, we require beyond the conjugation lemma the further asymptotic ODE tools 
of the convergence and tracking/reduction lemmas of [6S \ I44| . Statements and proofs of these 
results are given for completeness in Appendix [Xl 



4.1 The small-a and -fiQ limits 



Proof of Theorem \ 1.1b] (a ^ 0). Considering (|2.1|) as indexed hy p := a with A = we 
have ()A.2p - ()A.3P by uniform exponential convergence of £> as x ^ ±00. Take without loss of 
generality /i = 1. Applying Lemma lA.H we find that the transformations conjugating (12. ip 
to its constant-coefficient limits Z' = A'^Z, by which the Evans function is defined in (j2.18p . 
are given to 0{a) by the transformations P^ conjugating to its constant-coefficient limits the 
(7 = system W = A°(x, X)W, with 



A^ 
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\v 


V 




















Vo 


-Blv 


Xv 


0/ 



(4.1) 



that is, P^ = P^ + 0{a). Moreover, for bounded from zero, and cr/A sufficiently small, it is 
straightforward to verify that the stable subspace of A+ (A, a) is given to order cr/A by the span of 



and (0, q+Y' , where is the stable eigenvector of 



1 

\v V 



and q+ = {-^a^XQ/X, 1), 

T 



and, similarly, the unstable subspace of A^ (A, a) is given to order cr/A by the span of (r__ , s 
and (0, where is the unstable eigenvector of 1^ and = {y^aj^io/X, 1). Thus, 

the Evans function for a > 0, appropriately rescaled, is within 0{a/X) of the product of the 
Evans function of the diagonal block w' = ^ w initialized in the usual way, which is 

nonzero by our earlier analysis of the decoupled case = 0, and of the trivial flow 



w 



vixj 




A 



w 



(4.2) 



initialized with vectors parallel to q^^ in the conjugated flow. To estimate the second deter- 
minant, we produce explicit conjugators P^ for the a = flow, making use of the observation 
of Remark 1 2 . 1 01 that . by lower block-triangular form of the A^, these may be taken lower block- 
triangular as well, and so the problem reduces to finding conjugators for the flow (14. 2p in 
the lower block. But, these may be found by exponentiation to be 



1 

c±{x) 



(v - v±){y)dy, 
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yielding an Evans function to order a/A of 



det(p_g_,p+g+) - det ^ ,_(o)y^ i _ c+ ' 

or 

(1 + 0(v^))2 V(T//o/A / 0. 
In particular, the Evans function is nonvanishing for 1 S> \ y^a/X\ » |cj|, |<t/A|, as occurs for 

cr< |A| < a'^. 

Since the Evans function is nonvanishing in any case for |A| sufficiently large, by Theorem 11.9^ 
we obtain nonvanishing except in the case < |A| < Ccr, which must be treated separately. 

To treat |A| < Ca, notice that the stable/unstable subspaces of decouple to order a for 
cr > sufficiently small and A only bounded into the direct sum of (r, s)^ already discussed and 
(0)^)"^ with the stable/unstable eigenvectors of 

/ 1 A 

\X/afio a^QvJ ' 

which may be chosen holomorphically as q± = =F \/ct2/Xq/4 + X/afio), with a single 

square-root singularity at A = — (cT/io)^/4. For |A| > Ca, these factor as 

q± = (1 + 0{a/X)){l, TVWf^o) = (1 + 0{a/ X)){y^X/^o)q±, 

and the Evans function for ex > by our previous computations thus satisfies 

D'^iX) = (A/cT^o) X (1 + 0(y^))2e^-(°)+^+(°)\/^^ ~ C^/X/^o. 

Taking the winding number about |A| = Ca on the punctured Riemann surface obtained by 
circling twice the branch singularity A,, = (cr^o)^/4, we thus obtain winding number one. 
Subtracting the nonnegative winding number obtained by circling twice infinitesimally close to 
A*, we find (similarly as in the treatment of the large- amplitude limit, case Bf > yfjlo) that 
there is at most one root of on 5RA > 0, for u > sufficiently small, so that stability is decided 
by the sign of the stability index, which is the product of a transversality coefficient and the 
hyperbolic stability determinant (resp. low-frequency stability condition, in the overcompressive 
case). A singular perturbation analysis of (|B.4p as o" ^ shows that (since it decouples into 
scalar fibers) connections are always transverse for o" > 0, so the transversality coefficient does 
not vanish. Hyperbolic stability holds always for Lax 1- and 3-shocks (Proposition II .41) . and the 
low-frequency stability condition holds for intermediate (overcompressive) shocks by a similar 
singular perturbation analysis, so the stability determinant does not vanish either. 

Thus, the sign of the stability index is constant, and so there is always either a single 
unstable root of D" on 5RA > or none, in each of the three cases. But, the former possibility 
may be ruled out by homotopy to the stable, small-amplitude limiting case. Thus, all type 
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shocks are reduced Evans stable for a > sufficiently small. The asserted ^/a asymptotics 
follow from the estimates already obtained in the proof; uniform convergence to zero follows 
by estimating instead to order \/o", at which level we obtain a determinant involving two 
copies of the constant solution W = (0, 0, 0, 1)^ of the limiting £7 = system, giving zero as the 
limiting value. □ 



Proof of Theorem \1.16\ (^o 0). The case /io — > is similar to but a bit tricker than the case 
o" — > just discussed. Fixing without loss of generality = 1 and applying Lemma lA.ll 
we deduce that the transformations Pl^° conjugating (12. ip to its limiting constant-coefficient 
systems, by which the Evans function is defined in (j2.18p . satisfy Pl^" = P^ + O{fj,o), where P^ 
are the transformations conjugating to its constant-coefficient limits the upper block-triangular 
^0 = system 



(4.3) 



which has a constant right zero-eigenvector r = {aB^,0, X,0)^ and an orthogonal constant 
left zero-eigenvector £ = (0,0,0, 1)"^, signaling a Jordan block at eigenvalue zero. It is readily 
checked for the limiting matrices at ±00, similarly as in the o" — > case, that for ^o/X sufficiently 
small, the Jordan block splits to order ~ fio/X, so that the "slow" stable eigenvector at -|-oo 
(that is, the one with eigenvalue near zero) is given by 
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r + c+y^VA(*, *, *, 1)^ + O{fio/X), 
and the slow unstable eigenvector at —00 by 



r + c_ vWA(*, *, *, 1)^ + O(^o/A), 

where c± are constants with a common sign. (Here, we deduce nonvanishing of the final 
coordinate of the second summand without computation by noting that the dot product with 
i must be ~ ^/JIoJX.) 

As for the u ^ case, we now observe that (|4.3p may be conjugated to constant-coefficients 



by block-triangular conjugators P± 
system 



1 



, where p± conjugate the upper lefthand block 



1 

■0 

-B\v 

Moreover, changing coordinates to lower block-triangular form 





(4.4) 



Wi 

Mi - XWi/aBl^ 



1 

V 



-Bp - X/aBl 




(4.5) 
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conjugating by a lower block-triangular conjugator, and changing back to the original coordi- 
nates, we see that the conjugators p± may be chosen to preserve the exact solution {Wi, W2, W4, W^)'^ 
r. 

Combining these observations, we find that the Evans function for A bounded and Ho/X 
sufficiently small is given by 

A ,f 4 ^ ^ , /\A (4-6) 

= det , — — ^ ^ , — — + Oifio \) 

= (c_ + c+) vWAd(A) + 0(/io/A), 

where r =: and d{X) := det (wg^ f) is a nonstandard Evans function associated with 

the upper-block system (j4.4p . where and v'^ as usual are unstable and stable eigendirections 
of the coefficient matrix, but we have included also the neutral mode r. Expressed in coordinate 
(j4.5p . d{X) reduces, finally, to the standard Evans function d{X) for the reduced system 

Wa - XWi/aBl) ~ \-Blv - X/aB^ J - XWi/aB^ 
which may be rewritten second order equation 

{X + a{Blf)z + z' = iz'/vy (4.7) 
in 2; = Taking the real part of the complex L^-inner product of z against (|4.7p gives 

5f?A||^;||i. + \\v^{^X + aiBir)\\l, = -\\v'/Vd\\h, 

contradicting the existence of a decaying solution for SRA > and verifying that d{X) 7^ 0. 
Consulting (j4.6p . therefore, we find that D^^"{X) for A bounded and fio/X sufficiently small does 
not vanish and, moreover, ~ c^/JxqJX for A sufficiently small, c 7^ constant. Performing 
a Riemann surface winding number computation like that for the case o" — > 0, we find, finally, 
that does not vanish for any 3f?A > 0. We omit the details of this last step, since they 
are essentially identical to those in the previous case. Likewise, the asserted asymptotics follow 
exactly as before. □ 



4.2 The large- and small-_B* limits 

Proof of Theorem \1.15[ Stability in the small- -BJ' limit follows readily by continuity of the Evans 
function with respect to parameters, the high-frequency bound of Theorem 1 1.91 and the zero-B^ 
stability result of Proposition IB.ll We now turn to the large- -B| limit. Let us rearrange (j2.ip . 
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By Theorem 11.91 we have stabihty for |A| > C|i?Jp independent of For vj^ > 0, we 
find easily stability for |A| > Cji^J'l for B^ sufficiently large. For, rescaling x — > and 
W {X^/^Wi, X^/^W2, Ws, VF4)^, we obtain W' = Aw = AqW + 0{\BI\-^)W, where 



An 














X'/^v 



Ai/2 
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X^/'^afiov 
—av 
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(4. 



with X^/"^ := A^/^/i?*, and v = v{x/Bl), v independent of 

For A » l-B^I^^ it is readily calculated that Aq has spectral gap » for SRA > 0. 

Indeed, splitting into cases A > and A <C 1, it is readily verified in the first case by 

standard matrix perturbation theory that there exist matrices R{v{x)) and L = R"^, both 
smooth functions of v, such that 



LAqR = D:-- 



M 
N 



with > 6 > and < —6 < 0. Making the change of coordinates W 
the approximately block-diagonal equations Z' = AZ, where 



A 



LAR - L'R = D + 0{\BI\~^). 



RZ, we obtain 



(4.10) 



Using the tracking/reduction lemma. Lemma \AA\ we find that there exist analytic functions 
Z2 = ^2{zi) = 0{r) and zi = $1(^2) = 0{r) such that {zi,^2{zi) and ($1(22), -22) are invariant 
under the flow of (j4.16p . hence represent decoupled stable and unstable manifolds of the flow. 
But, this implies that the Evans function is nonvanishing on A G {KA > 0} for B^ sufficiently 
large and |A|^/^ > \Bl\/C, for any fixed C > 0. See [55l|38l[50] for similar arguments. 

If |A^/^| <C Bl on the other hand, or, equivalently, |A^/^| <^ 1, then we can decompose A 
alternatively as W' = AW = BqW + X^/^BiW + 0(|S||-i), where 
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\0 -V 
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(4.11) 



By smallness of X^^'^ together with spectral separation between the diagonal blocks of Bq, 
there exist L, R, LR = such that the transformation W = RZ takes the system to Z' = 
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{LBR - L'R)Z = CZ + 0{\BI\-^)Q, wheie e = *^ and 

Diagonalizing /J into growing and decaying mode by a further transformation, and applying the 
tracking lemma again, we may decouple the equations into a scalar uniformly-growing mode, a 
scalar uniformly-decaying mode, and a 2-dimensional mode governed by 

z' = -\p-^z + 0{\BI\-^ + \X'^\)z. (4.13) 

If A ^ l-BJl^^, or, equivalently, |A| ^ \Bl\, then we make a further transformation diagonalizing 
p-^ at the expense of an 0(|i3^| ) error, then use the resulting > |A| spectral gap together 
with the tracking lemma to again conclude nonvanishing of the Evans function. 

Thus, we may restrict to the case |A| < Cj-BJ"!, or |A| < C|i?5[|~^- Considering again ()4.13p 
in this case, we find that all 0{\B^\~'^ + |A^|) entries converge at rate 

0{\Bl\-^)\v - v+\ < C|5i*r2e-l^l/^^i' 

to limiting values, whence, by the convergence lemma. Lemma lA.ll the Evans function for the 
reduced system (j4.13p converges to that for z' = —XP~^z as B^ — > ooH But, this equation, 
written in original coordinates, is exactly the eigenvalue equation for the reduced inviscid system 



a J \-Bl J \a 

which may be shown stable by an energy estimate as in the case o" = 0. 

Finally, noting that the decoupled fast equations are independent of A to lowest order, we 
find for |A| bounded and B^ ^ oo that the Evans function (which decomposes into the product 
of the decoupled Evans functions) converges to a constant multiple of the Evans function for 
(j4.14p . For |A| < C, or A < C\Bl\~'^, however, we may apply to (j4.13p the convergence lemma. 
Lemma lA.H together with Remark IA.3| to see that the Evans function in fact converges to 
that for the piecewise constant-coefficient equations obtained by substituting for the coefficient 
matrix on x ^ its asymptotic values at ±00, that is, the determinant d := det(r^,r~), where 
r+ is the stable eigenvector of 

-Bl/fiov\ 
-Bl ) 

at +CXD and r~ is the unstable eigenvector at —00. Computing, we have = (1,=F\/AW^)^ 
where giving a constant limit d = ^fiov+ + as claimed. □ 



® Here, as in Remark lA.3l we are using the fact that also the stable/unstable eigenspaces at +00/— 00 converge 
to limits as |-Bi| — > 00. Together with convergence of the conjugators P^^ , this gives convergence of the Evans 
function by definition p.lSp . 
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4.3 The large-(j|Uo limit 



Proof of Theorem \1.17\ By Theorem II. 9| it is sufficient to treat the case |A| < CafiQ. Decom- 
pose (|2.ip . = 1, as W = RAq + Ai, where R := crj^iQ, X := X/R, and 
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Blv 0/ 



with \Aa\ < C. If |A| > 1/C > 0, then the lower 2x2 righthand block of has eigenvalues 

uniformly bounded from the eigenvalues zero of the upper lefthand 2x2 block. By 
standard matrix perturbation theory, therefore, there exist well-conditioned coordinate trans- 
formations L, R depending smoothly on v such that 



D := LAqR 
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Xv 


Vo 






v\/l 
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-vVll 



Making the coordinate transformation W = RZ, we obtain Z' = DZ + 0{1)Z. Applying 
the tracking lemma. Lemma IA.41 we reduce to a system of three decoupled equation, con- 
sisting of a uniformly growing scalar equation, a uniformly decaying scalar equation, and 

oj ^ ~^ 0{R^^)z. Rescaling by z := ^1/2^ ^' obtain 



a 2 X 2 equation z' 



1 
Xv 



RXv 0, 

+ 0{R^^/'^)y, which, by a second application of the tracking lemma, may 



be reduced to a pair of decoupled, uniformly growing/decaying scalar equations, thus com- 
pletely decoupling the original system into four growing/decaying scalar equations, from which 
we may conclude nonvanishing of the Evans function. 

It remains to treat the case |A| ^ 1. We decompose ()2.ip . = 1, in this case as W = 
RBq + Bi , where 
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with |A| ^ R. Defining T 



I 

I 



where 9 = (0, —Blv/fio,v)'^, and making the change of 



variables W = TZ, we obtain Z' = RCqZ + CiZ, where 

A 



Co 



'Co 





(fill! 

MO 



+ A' 



Co = /?■ 



1 

Xv v-{Bl)^vy^,o XiBDHyfio 
-B^v^ A{)2 
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and 

/O 1 0\ . . 

P={Xv V 0], x={0 -Blv \v), Ci = ( *)=0(1). 

Vo 0/ ^* 

Applying the tracking lemma, we reduce to a decoupled system consisting of a uniformly 
growing scalar equation 

y' = {R+^-^lL + x)y^y + 0{l/R)y (4.15) 
associated with the lower right diagonal entry and a 3 x 3 system 

z' = coz + 0{l/R)z. 
/ 1 * \ 

For |A| 1, we may write cq = I * I +0(1), and apply the tracking lemma again to 

V Xv^l 

obtain three decoupled equations uniformly growing/decaying at rates ±\/A^ and At)^, giving 
nonvanishing of the Evans function. For |A| < C on the other hand, we may apply the con- 
vergence lemma, Lemma lA.H using the fact that the 0(1/ R) coefficient converges to its limits 
as CR~^e~^^^^, r/ > 0, together with Remark I A. 31 to obtain convergence to the unperturbed 
system z' = cqz. But, this may be recognized as exactly the formal limiting system (jl.2ip 
for {a = oo), which is stable by Theorem 11.71 (established by energy estimates). Noting that 
the Evans function for the full system is the product of the Evans functions of its decoupled 
components, and that The Evans function for the scalar component converges likewise to that 
for y' = {R + ^^j^ — h X)v'^y, or (by direct computation/exponentiation) d{X) = e'^oR+ci+c2X £qj. 
constants cj, we find, finally, that the full Evans function after renormalization by factor e~'^°^ 
converges to a constant multiple of the Evans function for (jl.2ip . □ 



4.4 The limit as /i/ (2yU + r^) — > or — cx3 

Finally, we briefly discuss the effect of dropping the gas-dynamical assumption rj = —4fi/3, 
and considering more general values of (2/u + rj) > 0. This parameter does not appear in the 
transverse equations, so enters only indirectly to our analysis, through its effect on the gas- 
dynamical profile v{x). Specifically, denoting r := ///(2// + ry) — > 0, and taking as usual the 
normalization /i = 1, we find that 

v{x) = v(rx), 

where v is a profile independent of the value of r. Thus, the study in [13] of the limit r — > 
is the limit of slowly-varying coefficients, and the opposite limit r — > oo is the limit rapidly- 
varying coefficients. We consider each of these limiting cases in turn. Intermediate values of 
/i/(2/i+r/) would presumably need to be studied numerically, an interesting direction for further 
investigation. 
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4.4.1 The r ^ limit 

In the r — > Umit, we have the following result completing the analysis of [13] 

Proposition 4.1. Parallel isentropic MHD shocks with ideal gas equation of state are reduced 
Evans stable in the limit as r ^ with other parameters held fixed. 

Proof. The case Bl < including Lax 1-type, overcompressive type, and some Lax 3-type 

shocks has been established in [13J by energy estimates. Thus, it suffices to treat the case of 
Lax 3-shocks and (by Proposition II. 9p bounded |A|. 

For shocks of any type, it is straightforward to verify that the Evans function is nonvanishing 
on A E {5RA > 0} \ B{0,e), any e > 0, for r sufficiently small. For, on this set of A, there is a 
uniform spectral gap between the real parts of the stable and unstable eigenvalues of A{x, A), 
for all X £ (— oo,+oo), by the hyperbolic-parabolic structure of the equations, similarly as in 
Lemma [2.1[ It follows by standard matrix perturbation theory that there exist matrices R{v{x)) 
and L = such that 

'M ' 



LAR = D :-- 



N 



with > 6 > and < —6 < 0. Making the change of coordinates W = RZ, we obtain 
the approximately block-diagonal equations 

Z' = AZ, (4.16) 

where 

A := LAR -L'R = D + 0{v^) = D + 0{rv^). (4.17) 

Using the tracking/reduction lemma, Lemma lA.41 we find that there exist analytic functions 
Z2 = ^2{zi) = 0{r) and zi = $1(2:2) = 0{r) such that {zi,^2{zi) and {^i{z2), Z2) are invariant 
under the flow of (I4.16P , hence represent decoupled stable and unstable manifolds of the flow. 
But, this implies that the Evans function is nonvanishing on A G {5RA > 0} \ .6(0, e), any e > 0, 
for r sufficiently small. See \55\ [38l [50] for similar arguments. 

Now, restrict to the case of a Lax 3-shock for A G {5RA > 0} H -6(0, e) and e > sufficiently 
small. By examination of A{x,X) at A = in the Lax 3-shock case, we find that on B(0,e) it 
has one eigenvalue //+ that is uniformly negative, one eigenvalue that is uniformly positive, 
and two that are small. By standard matrix perturbation theory [38[ 150] . there exist matrices 

L = \ Lq \ , R= {R+ Ro R-) with LR = I and L'jRj = such that 



LAR{x,X) = 

where the crucial factor A in AMq is found by explicit computation/Taylor expansion [55 
[50] 141], and Mq = + 0{X), where /3 as in (|2.5p is the hyperbolic convection matrix 

P:=(_l. -''\/^°'). (4.18) 
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Moreover, R, L depend only on v, A. Making the change of coordinates Z := LW, we obtain 
Z' = B{x, X)Z, where 

B = LAR -L'R= I 0{v^) XMq 0{v^) 




0{Vx) /i- 

Applying the tracking/reduction lemma again, we reduce to three decoupled equations associ- 
ated with the three diagonal blocks. The two scalar equations associated with fi± are uniformly 
growing/decaying, so do not support nontrivial decaying solutions at both infinities. Thus, 
vanishing of the Evans function reduces to vanishing or nonvanishing on the central block 

w' = (AMo + 0{vl))w, 

w e Noting that \\vI\\li = 0{r) 0, we may apply the convergence lemma, Lemma lA.lt 
together with Remark \A.3\ to reduce finally to 

z' = XMqz, Mo := LqARo. (4.19) 

For |A| < r, we have |AMo - AMo(+oo)| < CAe"^'', hence 

IIAMo - AMo(+cx))||ii[o,+oo) = 0{X/r) ^ 

and we may apply the conjugation lemma to obtain that the Evans function for the reduced 
central system (|4.19p is given by (1 + 0(A/r)) det(r_, r_|_), where r_ is an unstable eigenvalue 
of Mo(— oo) and r_|_ is a stable eigenvalue of Mo(— oo). Noting that these to order A are 
stable/unstable eigenvectors of P±, we find by direct computation that the determinant does 
not vanish. Indeed, this is exactly the computation that the Lopatinski determinant does not 
vanish for 3-shocks. Thus, we may conclude that the Evans function does not vanish for |A| <^ r. 

Finally, we consider the remaining case r/C < \X\ < Cr, for C > large but fixed. In this 
case, we may for the same reason drop terms of order A^ in the expansion of AMq, to reduce 
by an application of the convergence lemma, Lemma lA.H and Remark \A.3\ to consideration of 
the explicit system z' = —XP~^z, which is exactly the inviscid system 



a J \—Bl 1 y V" 

But, this may be shown stable by an energy estimate as in the case a = 0. Thus, we conclude 
that the Evans function does not vanish either for |A| ~ r and 3ftA > 0, completing the proof. □ 

Remark 4.2. The Lax 1-shock case may be treated by a similar but much simpler argument, 
since growing and decaying modes decouple into fast and slow modes. The overcompressive case 
is nontrivial from this point of view, since v passes through characteristic points as x is varied. 
However, we conjecture that the argument could be carried out in this case by separating off 
the single uniformly fast mode and treating the resulting 3- dimensional system by an energy 
estimate like that in the a = or ^ case. 
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4.4.2 The r ^ oo limit 

The opposite hmit r ^ cxd is that of rapidly- varying coefficients, and is much simpler to carry 
out. By the change of coordinates x x/r, we reduce v{x) to a uniformly exponentially 
decaying function v{x), and the coefficient matrix A{x, A) to a function A = r~^A that decays 
to its limits as 



where 9 > Oq > 0. Applying the convergence lemma, Lemma lA.H together with Remark \A.3\ 
we obtain the following simple result. 

Proposition 4.3. In the limit r oo, the reduced Evans function converges uniformly on 
compact subsets of > to D^{X) = det{R^ , R^), where are matrices solving Kato's 
ODE, whose columns span the stable (resp. unstable) subspaces of A±. 

That is, determination of stability reduces to evaluation of a purely linear algebraic quantity 
whose vanishing may be studied without reference to the evolution of a variable-coefficient ODE. 
This can be seen in the original coordinates by the formal limit 



We examine stability of D numerically, as it does not appear to be readily accessible analyti- 
cally. 

5 Numerical Investigation 

In this section, we discuss our approach to Evans function computation, which is used to 
determine whether any unstable eigenvalues exist in our system, particularly in the intermediate 
parameter range left uncovered by our analytical results in Section 11.71 Our approach follows 
the polar-coordinate method developed in ^2]; see also [H |271 [291 ISSl HI]- Since the Evans 
function is analytic in the region of interest, we can numerically compute its winding number 
in the right-half plane around a large semicircle B{0,A) n {3f?A > 0} containing (jl.22p . thus 
enclosing all possible unstable roots. This allows us to systematically locate roots (and hence 
unstable eigenvalues) within. As a result, spectral stability can be determined, and in the case 
of instability, one can produce bifurcation diagrams to illustrate and observe its onset. This 
approach was first used by Evans and Feroe [12] and has been applied to various systems since; 
see for example pi [2| [HI 17]. 

5.1 Approximation of the profile 

Following [11[27], we can compute the traveling wave profile using one of MATLAB's boundary- 
value solvers bvp4c [U], bvp5c [36], or bvp6c [21j, which are adaptive Lobatto quadrature 
schemes and can be interchanged for our purposes. These calculations are performed on a 



\A{x, X)-A±\< Cr-^e-^l^'l for x^O, 
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finite computational domain [— with projective boundary conditions M±{U — U±) = 0. 
The values of approximate plus and minus spatial infinity L± are determined experimentally 
by the requirement that the absolute error \U{zizL±) — U±\ be within a prescribed tolerance, 
say TOL = 10~^; see \29\ Section 5.3.4] for a complete discussion. Throughout much of the 
computation, we used L± = ±20, but for some rather extreme values in our parameter range, 
we had to lengthen our interval to maintain good error bounds. 

5.2 Approximation of the Evans function 

Throughout our numerical study, we used the polar-coordinate method described in [32], which 
encodes W = where "angle" Q = loi A ■ ■ ■ A lv^ is the exterior product of an orthonormal 
basis {ujj} of Span{Wi, . . . ,Wk} evolving independently of p by some implementation (e.g., 
Drury's method) of continuous orthogonalization and "radius" p is a complex scalar evolving 
by a scalar ODE slaved to 0, related to Abel's formula for evolution of a full Wronskian; see [32] 
for further details. This might be called "analytic orthogonalization", as the main difference 
from standard continuous orthogonalization routines is that it restores the important property 
of analyticity of the Evans function by the introduction of the radial function p (0 by itself is 
not analytic); see [32^ I54j for a discussion on this method. 

5.2.1 Shooting and initialization 

The ODE calculations for individual A are carried out using MATLAB's ode45 routine, which 
is the adaptive 4th-order Runge-Kutta-Fehlberg method (RKF45). This method is known to 
have excellent accuracy with automatic error control. Typical runs involved roughly 300 mesh 
points per side, with error tolerance set to AbsTol = le-6 and RelTol = le-8. 

To produce analytically varying Evans function output, the initial data V(— L_) and V(L_|_) 
must be chosen analytically using (12. 6p . The algorithm of [9] works well for this purpose, as 
discussed further in [H |32] . 

5.2.2 Winding number computation 

We compute the winding number by varying values of A around the semicircle B{0, A)n{5RA > 0} 
along 120 points of the contour, with mesh size taken quadratic in modulus to concentrate 
sample points near the origin where angles change more quickly, and summing the resulting 
changes in arg(Z)(A)), using Qlog D{X) = argZ)(A)(mod27r), available in MATLAB by direct 
function calls. As a check on winding number accuracy, we test a posteriori that the change in 
argument of D for each step is less than 0.2, and add mesh points, as necessary to achieve this. 
Recall, by Rouche's Theorem, that accuracy is preserved so long as the argument varies by less 
than vr along each mesh interval. 

5.3 Description of experiments: broad range 

In our first numerical study, we covered a broad intermediate parameter range to demonstrate 
stability in the regions not amenable to our analytical results in Section 11.71 a-^d also to close 
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our study for unconditional stability for all (finite) system parameters. Since Evans function 
computation is essentially "embarrassingly parallel" , we were able to adapt our STABLAB code 
to take advantage of MATLAB's parallel computing toolbox, sending to each of 8 "workers" on 
our 8-core Power Macintosh workstation, different values of A producing a net speedup of over 
600%. The following parameter combinations were examined: 



{-f,v+,Bl,no,a) e {1.0,1.1,11/9,9/7,7/5,5/3,1.75,2.0,2.5,3.0} 

X {0.8, 0.6, 0.4, 0.2, 10"^ , 10"^^ 10~^ 10~^ 10"^} 
X {0.2,0.8,1.4,2.0,2.6,3.2,3.8} 
X {0.2,0.8,1.4,2.0,2.6,3.2,3.8} 
X {0.2,0.8,1.4,2.0,2.6,3.2,3.8}. 



In total, this is 30,870 contours, each consisting of at least 120 points in A. In all cases, we 
found the system to be Evans stable. Typical output is given in Figure [TJ 

We remark that the Evans function is symmetric under reflections along the real axis (con- 
jugation). Hence, we only needed to compute along half of the contour (usually 60 points in 
the first quadrant) to produce our results. 

5.4 Description of experiments: limiting parameters 

The purpose of our second study is to verify convergence in the large-amplitude limit (t>+ — > 0) , 
as well as illustrate the analytical results the limiting cases, namely as i?| — > oo, Bl 0, 
^0 ~^ cxD, fiQ ^ 0, a ^ oo, fj ^ 0, r ^ cxd, and r ^ 0. In all cases, we found our results to be 
consistent with stability. 

In Tabledl we provide typical relative errors between the normalized and limiting-normalized 
Evans functions in the large-amplitude limit; we varied Bl for illustrative purposes. The relative 
errors are given by computing, respectively, 



along the contours except for small A (that is, when |A| < 10 ^). Note that in the large- 
amplitude limit, the relative errors go to zero, as expected. 

A The convergence and tracking lemmas 

A.l The convergence lemma 

Consider a family of first-order equations 



max 

j 



and max 

j 



D{Xj)-D^{X,) 
D^iXj) 



W = AP{x, X)W 



(A.l) 
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B]^ = (J.z 


B\ = (J.o 


Bl = 1.4 


Bl = 2 


Bl = 2.6 


Bl = 3.2 


Bl = 3.8 


lO(-l) 


9.94(-l) 


1.23 


3.46 


9.33 


2.16(1) 


4.89(1) 


1.09(2) 


10(-2) 


4.36(-l) 


5.19(-1) 


1.36 


2.82 


4.92 


8.19 


1.32(1) 


10(-3) 


1.42(-1) 


1.72(-1) 


4.50(-l) 


8.34(-l) 


1.25 


1.86 


2.53 


10(-4) 


4.23(-2) 


5.04(-2) 


1.32(-1) 


2.30(-l) 


3.23(-l) 


4.55(-l) 


5.88(-l) 


10(-5) 


1.26(-2) 


1.50(-2) 


4.00(-2) 


6.83(-2) 


9.35(-2) 


1.28(-1) 


1.61(-1) 


10(-6) 


3.94(-3) 


4.77(-3) 


1.28(-2) 


2.18(-2) 


2.96(-2) 


4.03(-2) 


5.01(-2) 


10(-7) 


2.16(-3) 


2.62(-3) 


7.08(-3) 


1.20(-2) 


1.63(-2) 


2.21(-2) 


2.75(-2) 


10(-8) 


2.07(-3) 


2.51(-3) 


6.78(-3) 


1.15(-2) 


1.56(-2) 


2.12(-2) 


2.63(-2) 



Table 1: Relative errors for D{X) and D{X). Here a = fiQ = 0.8 and 7 = 5/3. 



indexed by a parameter p, and satisfying exponential convergence condition (j2.14p uniformly 
in p. Suppose further that 

\{AP - Al) - {A^ - Al)\ < Clple-^^''^ e>0 (A.2) 

and 

\{AP-A')^)\<C\p\. (A.3) 

Then, we have the following generalization of Lemma [2. 91 a simplified version of the convergence 
lemma of [44j . 

Lemma A.l. Assuming (j2.14p and ()A.2p - ()A.3p . for \p\ sufficiently small, there exist invertihle 
linear transformations P^{x, A) = / + Q^{x, A) and (x, A) = / + @^ (x, A) defined on x > 
and x <0, respectively, analytic in A as functions into L°°[0, ztcxo), such that 

\{PP - P°)±| < C7i|p|e-^>l for x ^ 0, (A.4) 

for any < 9 < 9, some Ci = Ci{9,9) > 0, and the change of coordinates W =: P^Z reduces 
(rO|) to 

Z' = AP{x, X)Z for x^O. (A.5) 

Proof. Applying the conjugating transformation W {P^)~^W for the p = equations, we 
may reduce to the case that A^ is constant, and = /, noting that the estimate (IA.2P persists 
under well-conditioned coordinate changes W = QZ, Q{zizoo) = I, transforming to 

I {Q-'APQ - Q-'Q' - AP^) - {Q-'A^Q - Q'^Q' - Al) \ 

< \Q{{AP - Al) - (^0 - Al))Q-^\ + \Q-'\AP - ^0)±Q - [A^ - ^0)±|, 

where 

\Q~\AP - A^)±Q - {AP - A0)±| = 0{\Q - I\)\{AP - A^)±\ = 0(e-^l^l)|p|. (A.7) 
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In this case, ()A.2p becomes just 

\AP-AP^\ < CiHe-^l^ 
and we obtain directly from the conjugation lemma, Lemma 12.91 the estimate 

\Pl - Pl\ = \Pl - I\ < CCi|p|e-^>l 
for X > 0, and similarly for x < 0, verifying the result]^ □ 

Remark A. 2. In the case = constant, or, equivalently, for which ()A.2p is replaced by 
\AP - < Cilple-^l"^, we find that the change of coordinates W = P^Z , P| := (P°)±^P|, 
converts (rO|) to Z' = A^Z, where P^ = I + with 

le^l < CCi|p|e-^>l. (A.8) 

That is, we may conjugate not only to constant- coefficient equations, but also to exponentially 
convergent variable- coefficient equations, with sharp rate (jA.Sp . 

Remark A. 3. As observed in provided that the stable /unstable subspaces of A^/A^ con- 
verge to those of A^/A^, as typically holds given ()A.3p - in particular, this holds by stan- 
dard matrix perturbation theory J34'l if the stable and unstable eigenvalues of A^ are spectrally 
separated" ()A.4p gives immediately convergence of the Evans functions to on compact 
sets of X, by definition (|2.18p . 

A. 2 The tracking lemma 

Consider an approximately block-diagonal system 

W = (^^Q^ {x,p) + 6{x,p)Q{x,p)W, (A.9) 

where is a uniformly bounded matrix, 5{x) scalar, and p a vector of parameters, satisfying a 
pointwise spectral gap condition 

mincr(3f?Mf ) - max o-(5RM|) > r/(x) for all x. (A. 10) 

(Here as usual := (1/2) (A^ + N*) denotes the "real", or symmetric part of N.) Then, we 
have the following tracking/reduction lemma of [38l 



Lemma A. 4 ([Ml Consider a system ()A.9P under the gap assumption (jA.lOp . with 

uniformly bounded and rj S If sup {5/ t]){x) is sufficiently small, then there exist (unique) 

linear transformations ^i{x,p) and #2(2;, p), possessing the same regularity with respect to p as 



^The inclusion of assumption (|A.3|l . needed in (|A.7|) . repairs a minor omission in [44]. (It is satisfied for the 
applications in [44] . but is not listed as a hypothesis.) 
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do coefficients Mj and 6@, for which the graphs {{Zi, <I>2^i)} o^c^ {{^1(22), Z2)} are invariant 
under the flow of (jA.Op . and satisfy 

sup|$i|, sup|<i>2| < Csup{6/ri) (A. 11) 

and 

M{x)\<C e^>^''^^'5{y)dy, \^l{x)\ < C T 5{y)dy. (A.12) 

J X J —00 

Proof. By the change of coordinates x^x,5^5:=5/7] with dx/dx = ri{x), we may reduce 
to the case rj = constant = 1 treated in [38]. Dropping tildes and setting <I>2 := '^2'4^i , where 
ii^i, il^lY satisfies (1X9]) . we find after a brief calculation that <I>2 satisfies 

$'2 = (M2$2 - ^2Mi) + 5Q($2), (A.13) 

where Q is the quadratic matrix polynomial Q{^) := 621 + ©22*1' — <I?©ii + <^0i2^*. Viewed as 
a vector equation, this has the form 

^'^ = M<i>2 + SQ{^2), (A.14) 

with linear operator := M2^ — $Mi. Note that a basis of solutions of the decoupled 

equation <!>' = may be obtained as the tensor product $ = <j)(p* of bases of solutions of 
(j)' = M2(p and (p' = —M^cj), whence we obtain from (jA.lOj) 

e^' < Ce-''^ for z > 0, (A.15) 

or uniform exponentially decay in the forward direction. 

Thus, assuming only that $2 is bounded at —00, we obtain by Duhamel's principle the 
integral fixed-point equation 

$2(x) = r<I>2(x) := r e^^^-y^5{y)Q{^2){y)dy. (A.16) 

J ~oo 

Using ()A.15p . we find that T is a contraction of order 0{6/r]), hence ()A.16P determines a 
unique solution for 5/r] sufficiently small, which, moreover, is order 5/r] as claimed. Finally, 
substituting Q($) = 0(1 + |$|) = 0(1) in ([AlB]) . we obtain 

\Mx)\<C j ei^^-y^5{y)dy 

in x coordinates, or, in the original x-coordinates, (|A.12p . A symmetric argument establishes 
existence of $1 with the asserted bounds. Regularity with respect to parameters is inherited as 
usual through the fixed-point construction via the Implicit Function Theorem. □ 
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Remark A. 5. For r/ constant and 6 decaying at exponential rate strictly slower that e~^^ as 
X +00, we find from ()A.12p that ^2ix) decays like 6/r] as x ^ +00, while if 6{x) merely 
decays monotonically as x ^ —00, we find that ^2ix) decays like {6/rf) as x ^ —00, and 
symmetrically for <I>i. This and ()A.12p is a slight addition to the statement of \38l \44h which 
did not include pointwise information. We will not need this observation here, but record it for 
general reference/completeness. 

Remark A. 6. A closer look at the proof of Lemma \A.4\ shows that, in the approximately block 
lower-triangular case, 502i not necessarily small, there exists a block-triangularizing transfor- 
mation ^2 = 0{sup\6/r]\) « 1, under the much less restrictive conditions 

sup flV^Kieiil + 19221)) < 1 and sup(|5/r?l |02i I) « 



sup \6/ri\ 

(We do not use this here, but remark it for general application.) 

B Miscellaneous energy estimates 
B.l Stability for fi* = 

Proposition B.l. Parallel ideal gas MHD shocks are stable for B* = {) provided that the 
associated gas-dynamical shock is stable. 

Proof. For B\ = 0, the eigenvalue equations become 

\u + u = iiu' Iv, 

(B 1) 

Xa + a' = {l/afiQv){a' /v)' , 

or 

Xvu + vu = fj.u" , 

Xva + va = {l/afj,Q){a /ij)' . 

Taking the real part of the complex L'^-inner product of u against the first equation and a 
against the second equation and summing gives 

3f?A( y i)(|np + lap) = - j {^\u'\'^ + {I / a ^xov)\a' \^) + j vM'^ + l«P) < 0, 

a contradiction for 3f?A > and n, a not identically zero. Thus, we obtain spectral stability in 
transverse fields (tt, B) for i?* = so long as the profile density is decreasing Vx < 0, as holds 
in particular for the ideal gas case, either isentropic or nonisentropic. Likewise, transversality 
and inviscid stability criteria are easily verified in this case by the further decoupling of u and 
B equations. Stability in the decoupled parallel fields {v,ui) is of course equivalent to stability 
of the corresponding gas-dynamical shock. □ 

Remark B.2. By continuity, we obtain from the above also stability for magnetic field Bl 
sufficiently small. Stability for small magnetic field was already observed in (2^ \19^ . by a 
similar continuity argument. 
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B.2 Stability for infinite jj,o 

Proof of Theorem \1.7\ case fiQ. For ^0 = 00, equations p.20|) become 

Xw + w' = uw" /v, 
Xa + a — Biw =0, 

hence the w equation decouples and is stable by the argument for = 0. Thus, w = for 
3?A > 0, and so the second equation reduces to a constant-coefficient equation Xa + a' = 0, and 
thus is stable. □ 

B.3 Trans versafity for large 

Proposition B.3. For > + max | y'^^f^, | ; and all 1 > v+ > 0, profiles (neces- 
sarily Lax 3 -shocks) are transverse. 

Proof. For fi = 1, the (transverse part of the) linearized traveling- wave ODE is 

Transversality is equivalent to nonexistence of a nontrivial solution of (|B.4p . Taking the real 
part of the complex inner product of 

fxo VV -Bl\ ( u 



^ ' l/a/ioy V-^i ^ ) \B 



against both sides of (lB.4p . noting that 



^B) ' \-Bl V J \BJ / J 2 
and estimating |^| < ^ (see [27j, Appendix A for similar estimates), we obtain 

((|),*(M-^)(|)),0 



where := ( |j ^ ] and 



M — I ^0 /^^o \ f Ho -Bl 



Bl V J \ 1/afioJ \-Bl V 

is positive definite for B^ > The first minor of (M — N) is equal to the first minor of M, 

so positive for B\ > y/JiQ. Thus, M — A > 0, giving a contradiction, if B* > yJJIo and 

< det(M -N) = a[fiov - {Blf)^ - (/xq + ct/.o(5i*)') ^ (B.5) 
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for all 1 > t) > f-i- > 0. Estimating 

[^^ov - {BlfY = {Bl - ^of{Bl + ^o? = {B*, - ^^?[{Blf + 

we find that (jB.SP holds for (B* - y/JIof > max{2^, yielding the result. □ 

Remark B.4. What makes this argument work is the strong separation as B^ ^ oo of grow- 
ing and decaying modes, as evidenced by strong hyperbolicity of the coefficient matrix on the 
righthand side of ()B.4p . It could be phrased alternatively in terms of the tracking lemma of 
Appendix \A.S[ Also related are the "transverse" estimates of \1 7[ [7 ^ . 

B.4 High-frequency bounds 

Proof of 1 1 . 91 Multiplying the first equation of (jl.20p by vw, integrating in x along M, and 
simplifying gives 

A / + / vw'iv + fi / = — - / aw. 



^0 JR 

Taking the real and imaginary parts, respectively, gives 



3ftA / v\w 

and 



/ vH'^-l I WxkP + M ! \w'\^ = '^^ [ aw (B.6) 



/ v\w\'^ + Q [ vw'w = — ^9 / aw. (B.7) 

JR JR fJ-O Jr 

Adding and simplifying, noting that JRz + |92;| < V^\z\ and Vx < 0, yields 

(5RA + |9A|) / v\w\'^ + n I \w'\^ < [ v\w'\\w\ + ^:^^ [ \a\\w\. 



Using Young's inequality, and noting that Vx ^ and w < 1, we have 

(3?A+|9A|) / v\w\' + ^i [ \w'\^< ^+£2^^) / 
Jr Jr \ /^o y Jr 



vlwl'^ 



(B.8) 



4ei Jr 4:£2fio 

Multiplying the second equation of (jl.20p by va, integrating in x along M, and simplifying gives 



A / ulap + / va'a + 



1 f la'P 



criJ-o Jr V 

Taking the real and imaginary parts, respectively, gives 
5RA / via 



Bl / vw'a. 
Jr 



[ vla]"^ -If Vx\a\'^ + — f = f vw'a (B.9) 
Jr. ^ Jr (^fJ-o Jr V Jr 
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and 



Jr Jr Jr 

Adding and simplifying, again noting that JRz + < \/2|-z| and Vx < 0, yields 

(5iA+|9A|) / v|ap + — / < / v\a\\a\+V2Bl [ v\w'\\a 
Using Young's inequality, and noting that u < 1, we have 

J 1 2 



(3f?A + |9A|) / v\af + — [ \^ < (es + V2£4Bl) [ 
Jr cr/^o Jr V ^ ^ Jr 



4e4 Jr 463 J^ v 

Adding C x (|RTT]1 to (ira yields 

(5RA + |9A|) / i)(|u;|2 + C|a|2)+/^ /" |u;'|2 + /" ^ 

Jr Jr Jm 



[ v\w\'^ + C (es + V2e4Bl) [ 
Jr ^ Jr 



< I £i + g2 ^ I / v\w\^ + C [63 + V2e4Bl] / v\a 



2 



1 V2BIC\ f , 2 I C \ f |a'i2 



By setting 



\ 4ei 464 / \ 4e2/Uo 4e3 



1 SJ'o- /ioo- -, BIC 



this becomes 

|2 , ^|„|2 



(5RA + |9A|) / v(|u;r + C|ar) 

ipt ptoC^ / 7r V 2 1^ J Jr 



2/x 

This inequality fails for all choices of t<;, a, whenever 



3f?A + |9A| > maxj— H , — - H i }. 



Setting 

yields the right-hand side of ()1.22p . 
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Corollary B.5. Any eigenvalue A of (|1.20p with nonnegative real part satisfies 

4 V /"/^o 

Proof. Adding (IB.6P to C x (lB.9p . noting that < 0, w < 1, and using Young's inequality 
yields 

C f la'P 



3fiA / v{\w\^ + C\a\^) + fi [ \w'\^ + — [ 



f ( ^ „'-,l,.,|2 , 1 W?\ , T3*r< f ( ^ „'-,l„|2 , 1 |„„/|2 



Setting 



< ^ / eiv\w\' + — '—^ + BIC / e2v\a\' + —\w' 
fJ-o Jr V iei V J 7r V 



£1 = ^ and £2 ~ 



(B.15) 



4C 4// 

together with (jB.lSh . yields the right-hand side of (lB.14p . □ 

C Kato basis near a branch point 



By straightforward computation, fJ.±{\) := =f(^/2 + s/rf/A + X and V± := (1,;U±(A))"^ are 
eigenvalues and eigenvectors of the matrix A of (|2.8p in Example 12.51 The associated Kato 
eigenvectors are determined uniquely, up to a constant factor independent of A, by the 
property that there exist corresponding left eigenvectors such that 

{V ■ V)"^ = constant, {V ■ V)^ = 0, (C.l) 

where " ' " denotes d/d\; see Lemma l2.3l fiii). 

Computing dual eigenvectors = (A + fi^)~^{X, n±) satisfying (V • V)^ = 1, and setting 
= c±V^, = V^/c-t, we find after a brief calculation that (jC.ip is equivalent to the 
complex ODE 

■V-V\±^ _ f fi 
■ V -v 

which may be solved by exponentiation, yielding the general solution 

c±{X) = C{ri'/A + \)-^'\ (C.3) 
Initializing at a fixed nonzero poin10 , without loss of generality c±{\) = 1, we obtain formula 



''in the numerics of Section O we typically initialize at A = 10. 
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